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Abstract 

We discuss the low-energy analysis of models involving quarks and four-fermion 
couplings. The relation with QCD and with other models of mesons and meson 
plus quarks at low energies is discussed. A short description of how the heat-kernel 
expansion can be used to get regularization independent information, is given. 

The anomaly within this class of models and a physical prescription to obtain 
the correct flavour anomaly while keeping as much of the VMD aspects as possible 
is discussed. The major part is the discussion within this framework of the order p 4 
action and of two and some three-point functions to all orders in momenta and quark 
masses. Some results on hadronic matrix elements are given. 
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1 Introduction 



The problems of dealing with the strong interaction at low and intermediate energies 
are well known. At short distance we can use perturbative Quantum Chromo Dynamics 
(QCD)|Q but due to asymptotic freedom this can no longer be done at low energies. The 
coupling constant there becomes too large. A general method, that is, however, extremely 
manpower and computer intensive, is using lattice gauge theory methods. An overview of 
this field can be found in the recent lectures by Sharpe[[| or in any of the proceedings of 
the annual lattice conferences. 

At very low energies we can use the methods of Chiral Perturbation Theory (CHPT). 
A good overview of the present state of the art here can be found in the DA$NE workshop 
report [Q. CHPT is a rigorous consequence of the symmetry pattern in QCD and its 
spontaneous breaking. Both perturbative QCD and CHPT are good theories in the sense 
that it is in principle possible to go to higher orders and calculate unambiguously. The 
size of the higher orders also gives an estimate of the expected accuracy of the result. A 
disadvantage of CHPT is that as soon as we start going beyond lowest order, the number 
of free parameters increases very rapidly, thus making calculations beyond the lowest few 
orders rather impractical. We would thus like to obtain these free parameters directly from 
QCD. 

This has so far been rather difficult to do. The reason is that all available approaches,like 
lattice QCD, QCD sum rules[|J, etc. , have problems with enforcing the correct chiral 
behaviour. We would also like to understand the physics behind the numbers from the 
lattice calculations in a more intuitive fashion. Therefore there is a need for some models 
that interpolate between QCD and CHPT. We will require that these models have the 
correct chiral symmetry behaviour. 

It should be kept in mind that these are models and not QCD. The hope is that these 
models will catch enough of the essential part of the behaviour of QCD at low energies 
that they can be useful. Two major classes exist, those with higher resonances than the 
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pseudoscalars included and staying at the hadronic level, or those with some kind of quarks. 
Both of these have their drawbacks. In the first case there still tends to be a large number of 
parameters and in the second case most models do not include confinement. Confinement 
is treated by explicitly looking at colour singlet observables only. The other drawback is 
inherent in the use of a model. It is not possible to systematically expand and get closer 
to the "true" answer. 

We will look at models including some kind of constituent quarks. The main motivation 
is that the standard constituent quark picture explains the hadron spectrum rather well. It 
has problems when interactions have to be included. It also tends to break chiral symmetry 
explicitly. Here we do not attempt to explain the hadron spectrum but instead focus on 
the few lowest lying states only. 

The class of models we will look at, is those where the fundamental Lagrangian contains 
quarks and sometimes also explicitly meson fields. There exists a whole set of these models 
of increasing sophistication. Models that are mainly for study of the spectrum like the bag 
model are not included. See || for a review of various aspects of this whole area. 

The lowest member of the hierarchy are the quark-loop models. Here the basic premise 
is that interactions of mesons proceed only via quark loops. The kinetic term for the mesons 
is added by hand. As a rule these models have some problems with chiral symmetry. In 
particular pointlike couplings of more than one meson to a quark-antiquark pair have to 
be added in order to be consistent. This goes under various names like bare-quark-loop 
model. A version that incorporates chiral symmetry correctly and also considers gluons 
is known as the Georgi-Manohar model 0. Another variation is to use the linear sigma 
model coupled to quarks. 

The next level is what I would call improved quark-loop models. Here also the kinetic 
terms of the mesons are generated by the quark loops. The degrees of freedom corre- 
sponding to the mesons still have to be added explicitly by hand. This leads to somewhat 
counterintuitive results when calculating loops of mesons 0]. This class started as integrat- 
ing the nonanomalous variation of the measure under axial transformations and its most 
recent member is known as the QCD effective action model ||, that reference also contains 
a rather exhaustive list of references to earlier work. 

The third level differs from the previous in that it starts with a Lagrangian which is 
purely fermionic and the hadronic fields are generated by the model itself. The simplest 
models here are those that add four-fermion interaction terms to the kinetic terms for the 
fermions. These are usually known as extended Nambu-Jona-Lasinio (ENJL)|| models. 
They have the advantage of being very economical in the total number of parameters and 
of generating the spontaneous breakdown of chiral symmetry by itself. The previous class 
of models has the latter put in by hand. Most of the remainder will be devoted to this 
class of models. A review of the more traditional way of treating this model can be found 
in 
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The most ambitious method has been to find a chirally symmetric solution to the 
Schwinger Dyson equations. These methods are typically plagued by instabilities in the 
solution of the equations. In the end they tend to be more or less like nonlocal ENJL 



models. They typically also have a lot of free parameters. A recent reference is [ITT 
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Some common features of all these models are that they contain a type of constituent 
quark mass and confinement is introduced by hand. The quarks are integrated out in 
favour of an effective action in terms of colourless fields only. The analysis also assumes 
keeping only the leading term in the expansion in the number of colours, l/iV c ||12||, only. 
This is not always explicitly stated but there are very few papers trying to go beyond the 
leading term. 

I will concentrate on the ENJL models since they are the simplest ones where the 
spontaneous symmetry breaking and the mesonic states are generated dynamically rather 
than put in by hand. Various arguments for this model in terms of QCD exist, see [|H| [13 



and section [|. A physics argument for the pointlike fermion interaction is that in lattice 
calculations the lowest glueball mass tends to be around 1 GeV. So correlations due to 
gluons below this scale might be suppressed. 

In this review no attempt was made to get a complete reference list. For this I refer 
to the more standard review [|It]] but let me give a few more background references. The 
original model was introduced as a simple dynamical model to understand the pions as 
Goldstone bosons from the spontaneously broken chiral symmetry as originally suggested 
by Nambu. After the advent of QCD there were various attempts at deriving such a model 
from QCD, see e.g. [0. Then the model lay dormant for some time till it was revived in 
the early eighties by Volkov, Ebert and collaborators Jl5|. At about the same time a more 
theoretical argument for these models was given in 
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A partial list of references where 
the phenomenological success of this model was shown is |14| to JO . 

There has also been some work on the NJL model on the lattice. This was mainly 
concerned with the attempt of finding a continuum limit (cut-off to infinity) . 

In the mean time a parallel development took place in the derivation of the Wess- 
Zumino-Witten term[53[] from quark models|P|]. This approach was then also used for 
the non-anomalous part of the effective action |46j to |49j. This can be found reviewed in 
45fl . This was later extended to include gluonic effects || and applied to nonleptonic weak 



matrix elements plj] , the higher order "anomalous" effective action |52[ and the ir + — 7r° 
mass difference 0. The requirement of propagating pseudoscalars that was found in the 
last reference provided an extra reason to go to purely fermionic models including chiral 
symmetry breaking and possibly gluonic effects. 

In [|13| the first step was taken by a low-energy expansion analysis of the extended 
Nambu- Jona-lasinio model. This was then extended to all orders in momenta for two- 
point functions in the chiral limit in [53] and with non-zero current quark masses in |54] . 
Some work along similar lines can be found in [55. 56 1 and [57) but without the emphasis 
on regularization independence. In |53 also the 7r + — 7r° mass difference was analyzed. 
The extensions to three point functions can be found in |54| and the extensions due to 
the anomaly were discussed in |58| and used for the 7r°7*7* vertex in [54]. This work 
was then extended to the Bk parameter |59|, |60| . A low-energy analysis of more vector and 
axial- vector meson processes was also performed |TT| and the application to the muon g — 2 
discussed In addition several talks about this work have been given p3] , |5"4] , p7| . 
It is this series of work that is reviewed in this Physics Reports. 

The report is organized as follows. In section ^] we discuss the Nambu- Jona-Lasinio 
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model, its various extensions and its connection with QCD. In Sect. |3] the occurrence of 
spontaneous chiral symmetry breaking is discussed. The next section, f|, is a short overview 
of the low-energy hadronic Lagrangians whose parameters we will try to understand in the 
context of the ENJL model. The relation of ENJL to other models is discussed in the next 
section. The regularization method used and the arguments behind the regularization 
independence of some of the results are given in section || Sect. |7| discusses the implemen- 
tation of the QCD anomalous Ward identities within this framework. This is essentially 



the discussion given in ||58|| . Then we reach the main results reviewed here. 

The low-energy expansion analysis is in Sect. [8], the extension to all orders in momenta 
and quark masses in the next section, while some three-point functions are discussed in 
Sect. pi]. Here there is also a more general discussion of the emergence of vector meson 
dominance (VMD) and a more general meson dominance in this class of models. 

Then we give a short overview of the results for nonleptonic matrix elements obtained so 
far. These calculations are among the most nontrivial uses of the ENJL model performed 



so far, Sect. |TT[ In the last section, we briefly recapitulate the main conclusions. The 
appendices contain the derivation of the ward identities at one-loop to all orders in momenta 
and masses and the explicit expressions for some of the one-loop functions needed. For a 



review of the heat kernel expansion I refer to J45| and to || q(J an d |[3j f° r the specific 
notation used in this report. 



2 The Nambu-Jona-Lasinio model and its possible con- 
nection with QCD 

In this section the arguments for the ENJL model as a low-energy approximation to the 
QCD Lagrangian are discussed. The different ways of looking at this model are also 
presented from a QCD viewpoint. 
The QCD Lagrangian is given by 

£°QCD = -1^ + ^(9, + ^)?. (1) 
4 a=l 

We restrict ourselves here to low energies so the quarks are the up, down and strange 
quarks, q = (q^oV and q — (u d s). The gluons in (|I|) are given by the gluon field matrix 
in the fundamental SU(N C = 3) co i our representation, 

N c~ l \(a) 

G^^T^W (2) 

a=l A 

with (?w the gluon field strength tensor 

Gjg = d,G^ - d v Gf - g s f abc GfG^\ (3) 
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and g s the colour coupling constant (a s = g^/Air). The fields v^, a M , s and p are 3x3 ma- 
trices in flavour space and denote respectively vector, axial-vector, scalar and pseudoscalar 
external fields. 

The coupling constant in this Lagrangian decreases with increasing energy scales. This 
is known as asymptotic freedom and is the reason why at short-distances we can use QCD 
perturbation theory. The other side of the coin is that at long distances the coupling 
constant becomes strong and leads to nonperturbative physics. This is generally known as 
infrared slavery and is probably also responsible for the phenomenon of confinement. The 
QCD Lagrangian and a review book can be found in Ref. [Q. 

The Lagrangian in (|l|) has a large classical symmetry. There is of course the gauge 
symmetry SU(N C ). In addition the different quark flavours are conserved leading to a 
U{\)y symmetry. The latter is increased to a U{3)y = U(l)y x SU{3)y symmetry if the 
three quark masses become equal. In addition for a zero quark mass there is an additional 
U(1)a for that flavour since for zero quark mass the Lagrangian does not couple the left and 
right handed combinations. For the case of zero quark masses the full classical symmetry 
of the lagrangian becomes 

SU(N c ) local x (SU(3) L x SU(3) R x U(l) v x U(l) A ) global . (4) 

Not all of these symmetries survive quantization. The U(1)a is explicitly broken by quan- 
tum effects^. This effect is known as the anomaly. 

The introduction of the external fields and a M allows for the global symmetries to 
be made local. The explicit transformations of the different fields are: 

Ql -> 9L{x)q L and q R -> g R (x)q R , 
r^^v^ + a^ -> g R r^g R + ig R d^g R , (5) 

and 

s + ip -> g R (s + ip)g[ . (6) 

Here gi,g R G SU(3)l x SU(3) r . The U(1)a is not a full symmetry. The quark masses are 
included in the Lagrangian via the scalar external field s 

m u 

external • (7) 

m s 

The number of colours, which is equal to three in the physical world, we have left free in 



order to use it as an expansion parameter |L2|]. We work in an expansion in inverse powers 



of the number of colours, 1/N C , where at the same time the QCD coupling constant is 

1 The SU(3)l x SU(3)r is also broken by the anomaly. These breaking effects are, however, not directly 
coupled to the strong interaction so they do not prevent the use of these symmetries in the same way as 
happens for the U(1)a- 
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scaled so that N c acs remains constant in the large N c limit. For a review of methods used 
to prove things in QCD, see fl~2[ . 



At low energies the coupling constant becomes strong and we cannot simply do a 
perturbation series in the above Lagrangian. The objects we will use to obtain physical 
observables is the generating functional of Green's functions of the vector, axial-vector, 
scalar and pseudoscalar external fields, v, a, s,p: 



e 



iF(v,a,s,p) _ /q 



Texp (i J d 4 xC QC D(x) 



. (8) 



The generating functional can be calculated in the path-integral formalism as: 

e ir(w) = \\ VG.VqDq exp [i J d 4 xC QCD (q, q, G; v, a, s,p)) 

= ~Z I VG » eXp ( _i / d4x \ G< $ G{a) ^) VqDqexv (i J d^xqiDq^j , (9) 

where D denotes the Dirac operator 

D = 7^ + iG») - i^( Vfl + 75 a M ) + i{s - t^p). (10) 

This generating functional is sufficient to be known at zero quark mass, since once it is 
known for zero quark mass, the identification of Eq. (|7|) makes sure it is known for nonzero 
quark masses. So if we know it as an expansion in external fields, we also know it as an 
expansion in quark masses and external fields. This is the basic premise underlying the 
formulation of Chiral Perturbation Theory of Gasser and Leutwyler|68|. 



The main assumption underlying the approach described in this report is to write the 
generating functional of Eq. (||) in a different way. At very low energies this can be done 
using Chiral Perturbation Theory: 

e ir(l, ' a ' s ' p) = VUexp (i J c^chpt) • (11) 

For an explanation of the symbols I refer to Sect. |j. 

This form of the generating functional can be used at low energies at the price of in- 
troducing a relatively large number of free parameters. We would therefore like to find 
an alternative way that can also be applied at low to intermediate energies and has fewer 
parameters. At present this involves making more assumptions about the low- and in- 
termediate energy behaviour of QCD than is inherent in using ([11]). One main approach 
is essentially to rewrite the generating functional in a functional integral form where the 
underlying degrees of freedom are still the quarks. There are various variations on this 
approach but we will replace (§) by: 



e iT(v,a,s,p) 



J "VG'lVqDq exp {% J d 4 x£ E NJL(g, q, G; v , a, s,p) \ . (12) 
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Here the integral over the gluonic degrees of freedom is either absent or only over low-energy 
gluons, see below. The Lagrangian in ([12]) is given by 



£ 



QCD 



r A x , rS,P , M Q 




with 



and 



£ 



S,P 
NJL 



N c Al 



a.b 



£ 



V,A 
NJL 



8'K 2 G V (A > 

N c Al 



Ytli&^qiWLWl) + (L ^ R) 

a,b 



(13) 



(14) 



(15) 



The couplings Gs and Gy are dimensionless quantities. For later convenience we also 
introduce the abbreviations 



9s 



4n 2 G s (A, 
N c Al 



and 



9v 



87T 2 G'y(A, 

iV c A2 



(16) 



Notice that in sections ^ and |8| the symbol gy is also used for the vector-two-pseudoscalar 
coupling. In principle the extra couplings in (|T3|) should be calculable in QCD as a function 
of A x and the QCD couplings. In practice this requires knowledge of the nonperturbative 
domain of QCD and we will determine all of the new parameters involved empirically. 

In the mean time it might be useful to see how this type of interaction could originate in 
QCD. This is illustrated in Fig. |. In Fig. |l]a the one-gluon-exchange interaction between 




Figure 1: (a) Conventional one-gluon exchange between two quark vertices in QCD. (b) 
Local effective four-quark interaction emerging from (a) with the replacement in eq. (|P7D. 
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quarks is shown. If we replace the propagator by it's short-distance part only we obtain a 
point like interaction. This is done via the regulator replacement, 



1 rVA 



Q 



2 



dre- rQ , (17) 



in the gluon propagator. This leads in the leading 1/N C limit to terms of the form (|14D 
and flTSP with the constraint 

G S = 4G V = ^N C . (18) 

7T 

This perturbative estimate of the extra couplings is of course only valid at short distances 
where perturbative QCD can be applied. A reliable calculation from QCD would require 
knowledge of all the higher orders. In particular, the anomalous dimensions of the two 
operators fll4[) and fll5|) are different so QCD can lead to different predictions for these 
operators already at the leading order in 1/N C . The constraint G$ = 4Gy has also other 
possible origins. In particular if we want to understand SU(6) in the baryon sector there 
should be spin independence of the constituent quark couplings. This leads precisely to 
this constraint. The couplings Gs and Gy are 0(1) in the large-iV c limit. This can also 
be seen in Eq. flI8|) . 

In general we could think of the Lagrangian of ENJL fllTf ) as being rooted in QCD 
by taking ([|) and performing the integral over gluons. The resulting effective action can 
then be expanded in terms of local operators of quark fields. Stopping at dimension 6 and 
leading order in the number of colours the Lagrangian is then precisely of the form ( |T~3j) 
but without any gluonic degrees of freedom. This is the standard picture of the ENJL 
model. An alternative view is that we integrate out the short-distance part of gluons and 
quarks and again expand the resulting effective action in local operators leaving only the 
leading terms in N c and dimensions. This again leads to a Lagrangian of the type ( |13D but 
this time with low-energy gluons. Several ways of looking at these gluons are possible but 
they are certainly not treatable as perturbative gluons. We will treat them as a way to 
describe the gluonic effects on the vacuum, i.e., we only keep their effects via the vacuum 
expectation values of gluonic operators. This is the point of view as taken in Ref . || . One 
of the results of the work reviewed here is that in the end the effects due to this gluonic 
vacuum expectation values are surprisingly small. 

Some alternative arguments on the basis of renormalons and QCD sum rules also 
exits||7D[. These arguments lead to the constraint 



G v = . (19) 
There is fact some work done on extensions of the Nambu-Jona-Lasinio model including 



higher order terms. Examples are the nonlocal NJL-models|7T|, Schwinger-Dyson type 



approaches [11] and models with some explicit higher order terms[72, 73|. These are terms 
that are suppressed by higher powers of 1/A X . 

Here we will keep only the first terms in order to keep the number of parameters down 
to a reasonable level. 
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It should be emphasized that this model does not include confinement. We will circum- 
vent this problem by only looking at observables that are explicitly colour singlets. The 
intermediate lines can in principle go on-shell above a certain energy. Mostly we will avoid 
this problem by working in the domain of Euclidean momenta and then doing an extrap- 
olation to the Minkowski domain using Chiral Perturbation Theory. The latter method is 
especially important in the treatment of nonleptonic decays in Sect. [11]. 

We work at the leading order in 1/N C throughout. At this order, as remarked above, 
the effects of U(1)a breaking due to the anomaly are absent. The other effects of the 
anomaly are still present like the two-photon decay of the 7r°. The underlying cause of 
this difference is that the strong coupling constant as also goes to zero in the large N c 
limit while the electromagnetic coupling does not. One effect of this limit is that nonet 
symmetry becomes exact, i.e., there is also a light pseudoscalar in the singlet channel or 
the 7/ is also light. Some discussions about effective lagrangians including the anomalous 
effect of U(1)a breaking can be found in[[73[]. A way of treating it in the context of the 
ENJL-model has been reviewed in Jl(J . 



The presence of the extra pointlike interactions in ( |T3[ ) has in fact some interesting 
consequences for the anomalous sector j58|. This is described in Sect. [7|. 

One more remark is needed here. We always implicitly assume that the quarks in (§) 
and (|T2D are identical. I.e. there are no other couplings of the external fields t> M , a M , s and 
p present. In the nonlocal models the presence of extra terms is already required by the 
chiral symmetry. This assumption should also be kept in mind when judging the results 
from the ENJL model. 



3 Spontaneous Chiral Symmetry Breaking in the NJL 
model 

The original paper of Nambu and Jona-Lasinio|| was in fact written to show the pion 
as a Goldstone boson and to provide an explicit model of spontaneous chiral symmetry 
breaking. All evidence point towards a spontaneous breaking of the axial symmetry by 
quark vacuum expectation values, (qq), in QCD. In the large iV c limit there exists a proof 
of this by Coleman and Witten|75|. Lattice gauge theory also finds agreement with this 



scenario |76| and a recent reevaluation of (qq) in Finite Energy Sum Rules [77] also gave a 
value consistent with the standard scenario. 

In the Nambu- Jona-Lasinio model we first have to calculate the fermion propagator to 
leading order in 1/N C . This can be done via the Schwinger-Dyson resummation of graphs 
depicted in Fig. |2]. There is no wave function renormalization to this order in 1/N C and the 
mass can be self-consistently determined from the Schwinger-Dyson equation. This leads 
to the condition 

Mi = 7^-^(01:^:10), (20) 
d 4 p i 



■ Mi ■ |0> = (m) = (qq)i = -NAM, J 



(2vr) 4 p 2 - M? 
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Figure 2: Schwinger-Dyson equation for the quark propagator, which leads to the gap 
equation in eq. (p0|). 
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Figure 3: Plot of the dependence of the constituent quark mass Mj as a function of Gs for 
several values of rrii 



N c 



16tt z 



AMf r (-1, e,) 



(21) 



Also to this order in N c , the constituent quark mass of flavour i, Mi is independent of 
the momenta and only a function of Gs, A x and rrii, the current mass of the zth flavour 
quark. The dependence on Gs is via gs defined in ([16]). It is not dependent on Gy. The 
T function in Eq. fl2"T|) is a consequence of our regularization scheme (see Sect. |B|). 

The scalar quark- ant iquark one-point function (quark condensate) obtains a non-trivial 
nonzero value. This nonzero value breaks chiral symmetry spontaneously leading to the 
occurrence of a nonet of pseudoscalar Goldstone bosons. 

The dependence on the current quark-mass is somewhat obscured in eq. (pl|). The 
quantity ti appearing in ( pTj) is Mf/A^. In figure [|we have plotted the dependence of Mj 
on Gs for various values of m,- and A v = 1.160 GeV. It can be seen that the value of M, 
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for small rrii converges smoothly towards the value in the chiral limit for the spontaneously 
broken phase. This is an indication that an expansion in the quark masses as Chiral 
Perturbation Theory assumes for QCD is also valid in this model. However, it can also 
be seen that the validity of this expansion breaks down quickly and for rrii ~ 200 MeV 
we already have 2Mj ~ A x . We note that the ratio of vacuum expectation values for light 
quark flavours increases with increasing current quark mass at p 2 = in this model and 
starts to saturate for rrii > 200 MeV. In standard x?T this ratio is taken to be 1 at lowest 
order and its behaviour with the current quark mass is governed (at 0(p A )) by the following 
combination of coupling constants 2L$ + H2 [^| in the large N c limit. Expanding (pl|) in 
powers of rrii thus gives a prediction for this combination of parameters, see Sect. |8|. 

In effect, the inclusion of gluonic corrections for this case is known to order (G 3 ), see 



Ref. 13 



4 Low Energy Hadronic Lagrangians 



As discussed in Sect. |3| the SU(3)l x SU(3)r symmetry in flavour space is expected to 
be spontaneously broken down to SU(3)v in QCD. According to Goldstone's theorem, 
there appears then an octet of massless pseudoscalar particles (tt, K, rf). The fields of these 
particles can be conveniently collected in a 3 x 3 unitary matrix U(<&) with det£7 = 1. 
Under local chiral transformations 

U(x) - g R U(x)gl (22) 
Whenever necessary, a useful parametrization for U(<&), which we shall adopt, is 

-$(x)\ 



U (<&) = exp -iV2- 



fo 



(23) 



where f ~f n = 93.2 MeV and (A are Gell-Mann's SU(3) matrices with trA a Aft = 2<5 ab ) 



4-3 

V2 



( 



V2 



+ 



V 

V6 



71 



71 



— 7T 

V2 



+ -3- 



K 



K + 
K° 
-2-2- 



(24) 



The 0~ octet <f>(x) is the ground state of the QCD hadronic spectrum. There is a mass 
gap from the ground state to the first massive multiplets with 1~, 1 + and + quantum 
numbers. The basic idea of the effective chiral Lagrangian approach is that, in order to 
describe physics of the strong interactions at low energies, it may prove more convenient to 
replace QCD by an effective field theory which directly involves the pseudoscalar - octet 
fields; and, perhaps, the fields of the first massive multiplets 1~, 1 + and + as well. Since 
we work here in the leading order in 1/N C we have to add the singlet components as well. 
In particular we have to add to $ 

/ rf 

rf I . (25) 
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The chiral symmetry of the underlying QCD theory implies that V(v, a, s,p) in eq. (§) 
admits a low-energy representation 



e iT(v,a,s,p) 



— [ VUVSVV.VA,,^ f d4xC *ff {u ' s ' v ^' v ' a ' s ' p) 
z J M " 

— [ DUe i f d4 ' xCeff ( u ' v ' a ' s ' p ) , (26) 
Z J 



Z 

where the fields S(x), V^(x) and A^(x) are those associated with the lowest massive scalar, 
vector and axial- vector particle states of the hadronic spectrum. Both C^ff and £ e // are 
local Lagrangians, which contain in principle an infinite number of terms. The hope is that, 
for energies sufficiently small with respect to the spontaneous chiral symmetry breaking 
scale A x , the restriction of £^ and/or C e ff to a few terms with the lowest chiral dimension 
should provide a sufficiently accurate description of the low-energy physics. The success 
of this approach at the phenomenological level is by now confirmed by many examplesQ 
We will later derive the effective Lagrangians and C e ff from the Nambu-Jona-Lasinio 
cut-off version of QCD. 

Let us now briefly summarize what is known at present about the low-energy mesonic 
Lagrangians and £ e // from the chiral invariance properties of Cqcd alone. 

The terms in C e ff with the lowest chiral dimension, i.e. 0(p 2 ), are 

4?/ = \ti{trD,UD^ + trfrrf + rf x )} (27) 
where denotes the covariant derivative 

D^U = 8^ -t{v^ + a^)U + tU{v^-a^) (28) 

X = 2B (s(x)+tp(x)). (29) 

The constants fo and Bq are not fixed by chiral symmetry requirements. The constant fo 
can be obtained from it — ► \w decay, and it is the same which appears in the normalization 



of the pseudoscalar field matrix U(Q) in fl24|), i.e. 



/„-/* = 93.2 MeV. (30) 
The constant Bq is related to the vacuum expectation value 

{m^)w=u4,s = -/o^o(l + 0{M)). (31) 
The terms in C e ff of 0(p 4 ) are also known. They have been classified by Gasser and 



2 For recent reviews see e.g. Refs. Q to JB1[. 
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Leutwyler fl69f:| 



4//^) = L^tiD^UW^U) 2 + L 2 ti(D^D u U)ti(D^D v U) 

+L£x{p il XJ^D»XJD v XJ^D v U) + L A ti(D^D»U)ti(x ] U + W X ) 

+Lsti[D^D»U(x ] U + W X )} + L 6 [tx( X ^U + W X )} 2 

+L 7 [tT( X U - Wx)? + LMX ] U X ] U + x^X^) 

-iL 9 tr(K U D^UD u W + FjfDpWDvU) + L w tr(W F^UF L ^) 

+HMK V Fr,u + Fl v F Lllu ) + H 2 trtfx), 
where Fl^v and Fr^ are the external field-strength tensors 

F Rlxu = dfj,r„ - duTfj, -i[r^ r v ] 
associated with the external left (Z M ) and right (r M ) field sources 



(32) 



(33) 



(34) 



The constants Li and Hi are again not fixed by chiral symmetry requirements. The Lj's 
were phenomenologically determined in Ref. [^|. Since then, £1,2,3 have been fixed more 
accurately using data from K14 1 82| . The phenomenological values of the Lj's that will 
be relevant for a comparison with our calculations, at a renormalization scale fi = M p = 
770 MeV, are collected in the first column of Table |l|. 

By contrast with C e f f , which only has pseudoscalar fields as physical degrees of freedom, 
the Lagrangian £^ involves chiral couplings of fields of massive 1~, 1 + and + states to 
the Goldstone fields. The general method to construct these couplings was described a long 
time ago in Ref. [83]. An explicit construction of the couplings for l - , 1 + and + fields 
can be found in Ref. P4| . As discussed in Ref. [B5], the choice of fields to describe chiral 



invariant couplings involving spin-1 particles is not unique and, when the vector modes are 
integrated out, leads to ambiguities in the context of chiral perturbation theory to 0{p A ) 
and higher. As shown in [S5||, these ambiguities are, however, removed when consistency 
with the short- distance behaviour of QCD is incorporated. The effective Lagrangian which 
we shall choose here to describe vector couplings corresponds to the so-called model II in 



Ref. [85 



In the NJL model it is of course obvious that the different representations for the meson 
fields should be identical since the original model is formulated in terms of fermions only. 
The choice of fields for the mesons is purely a matter of choice during the calculation. 



3 There are more terms in principle because of the presence of the singlet component as well. These 
have all zero coefficients at the leading order in 1/N C . 



14 



The wanted ingredient for a non-linear representation of the chiral SU(3)l X SU(3)r = 
G group when dealing with matter fields is the compensating SU(3)v transformation 
h(Q,gL,R) which appears under the action of the chiral group G on the coset represen- 
tative of the G/SU{3)v manifold, i.e. 

-> fe£(<J>)^($, = , 9L >R )Z{®)gl , (35) 

where £($)£($) = U in the chosen gauge. This defines the 3x3 matrix representation of 
the induced SU(3)y transformation. Denoting the various matter SU(3)y multiplets by 
R (octet) and Ri (singlets), the non-linear realization of G is given by 

R^h($,g LtR )Rh\$,g L , R ) 

R x -> R u (36) 

with the usual matrix notation for the octet 

1 8 

fl=^VA (t ¥ ) . (37) 

The vector field matrix V^(x) representing the SU(3)y octet of 1~ particles; the axial- 
vector field matrix A^{x) representing SU(3)y octet of 1 + particles; and the scalar field 
matrix S(x) representing SU(3)y octet of + particles are chosen to transform like R in 
eq. (H), i.e. (h = h($, g L , R )): 

V» -> hV^ ; -> ; 5 -> /iSTi*. (38) 

The procedure to construct now the lowest-order chiral Lagrangian C^j is to write 
down all possible invariant couplings to first non-trivial order in the chiral expansion, 
which are linear in the R fields, and to add of course the corresponding invariant kinetic 
couplings. It is convenient for this purpose to first set the list of possible tensor structures 
involving the R fields, which transform like R in eq. (|36| ) under the action of the chiral 
group G. Since the non-linear realization of G on the octet field R is local, one is led to 
define a covariant derivative 

d^R=d^R+[Y^R] , (39) 

with a connection 

r„ = \{$% - ifa + + t[dp - •(«„ - rf} (40) 
ensuring the transformation property 

d^R -> hd^Rh ] . (41) 
We can then define vector and axial-vector field strength tensors 

Vpv = - dvV^ and = d^A u - d u A^ , (42) 
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which also transform like R, i.e. 



V»v -> hV^tf and A„ v -> M^/il (43) 

There is a complementary list of terms that can be constructed with the coset represen- 
tative and which transform homogeneously, i.e. like R in fl36|) . If we restrict ourselves 
to terms of 0(p 2 ) at most, here is the list: 

^ = - i{v, + a„)]£ - Z[d, - i{v, - a M )]e f } = i^t^ = & (44) 

£,[i£,v and (45) 
X± = ^ ± (46) 

$ = ^ t ±$ t M- (47) 

Notice that in fl4*0|) does not transform homogeneously, but rather like an SU(3)v 
Yang-Mills field, i.e. 

r M -> w>t + /^t. (48) 

The most general Lagrangian to lowest non-trivial order in the chiral expansion 
is then obtained by adding to C^L in eq. ( f2~7|) the scalar Lagrangian 

C s = hi (d,Sd»S - M 2 S S 2 ) + c m tr (S X + ) + Qtr (S^) 5 (49) 
the vector Lagrangian 



C v = -Itr - 2M y V M V") - -L [/ytr ( V,J< + ^) + W tr (^f, CD 



and the axial- vector Lagrangian 



+ ••• , 

(50) 



C A = -^tr (A, U A^ - 2M\A^) - ^=f A tr (A,J^ U ) + • • • . (51) 

The dots in C v and C A stand for other 0(p 3 ) couplings which involve the vector field 
and axial- vector field N 1 instead of the field-strength tensors and A^ u . They have been 
classified in Ref. fl8~5| . As discussed there, they play no role in the determination of the 
0(p A ) Li couplings when the vector and axial- vector fields are integrated out. 

The masses My, Ms and M4 and the coupling constants c m , Cd, /y, 9v and /a are 
not fixed by chiral symmetry requirements. They can be determined phenomenologically, 
as was done in Ref. [53]. Since later on we shall calculate masses and couplings only 
in the chiral limit, we identify My, Ms and M4 to those of non-strange particles of the 
corresponding multiplets, i.e. 

M V = M P = 770 MeV ; M s = M ao = 983 MeV (52) 
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and 



M A = M ai = 1260 ± 30 MeV. 



The couplings fy and gy can then be determined from the decays p° 
respectively, with the result 



(53) 

e + e~ and p — ► ttti 



0.20 



and 



\g v \ = 0.090. 



(54) 



The decay a\ — > 7T7 fixes the coupling to 



0.097 ±0.022, 



(55) 



where the error is due to the experimental error in the determination of the partial width, 
r(ai — > 7T7) = (640 ±246) keV. For the scalar couplings c m and Cd, the decay rate ao — > tjtt 
only fixes the linear combination [84 



2m* 

C d ± 7i if) o o Cm 

M a - m v ~ m l 



(34.3 ± 3.3) MeV. 



(56) 



In confronting these results with theoretical predictions, one should keep in mind that they 
have not been corrected for the effects of chiral loop contributions. 

In addition in order to describe vector interactions beyond those that can be described 
by the above terms there are more terms possible. These do however not contribute to 
CHPT coefficients of order p 4 when integrated out. For a list of these terms see Ref. |JT] . 

I will now shortly review the different ways vector mesons tend to be implemented. A 
review can be found in [KB" 



There is the way of gauging the U(3)l X U(3)r symmetry by a set of vector meson 
fields, L M and These can be given a mass term without breaking the local symmetry 
by introducing the external fields Z M and r M defined above. To the Yang-Mills Lagrangian 
and the lowest order Lagrangian for the pseudoscalar mesons, with and in the 
covariant derivative now, we add a term of the form 



-m tr 



+ Ru 



(57) 



The mass m corresponds to the vector meson mass in the chiral limit and the axial-vector 
mass becomes different due to a partial Higgs mechanism, the field L M — mixes with 
the pseudoscalars. Including vector mesons only in this formalism requires sending the 
"bare" pion decay constant to infinity. This is often referred to as the gauged Yang-Mills 
formulation. 

A variation on the Yang-Mills principle is the hidden gauge formalism |87j . This formal- 
ism also allows for only the vector mesons to be included. There are more free parameters 
here than in the previous formalism at first sight but if one allows for higher order terms 



in both formalisms they are fully identical. This was proven in [57|. Removing the axial 
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vector mesons from the simplest gauged Yang-Mills version leads to the hidden gauge ver- 
sion (vectors only) with the extra constant a = 1. The usual VMD requirement has a = 2. 
This version, a = 1, also corresponds to Weinbergs original formulation of an effective 
Lagrangian for vectors and pions[|~~ 



One can also include the vector mesons in the general form as described by Callan et 
al. [jS3[] . This is the formulation described earlier in this section. There is also a version 



possible where and transform linearly under the chiral symmetry. This is the version 
that the ENJL model ends up with most simply. 

The last version is to use antisymmetric tensor fields to describe the (axial-)vector 



mesons. This was the formulation chosen in [34, 68 . This can be related to the other 



approaches by choosing the field strength rather than the bare field as the interpolating 
fields for the vectors. 

All of these formalisms can lead to identical physics by introducing extra pointlike pion 
couplings and higher order couplings as well. As such it is a matter of taste which version 
one chooses. Some of them tends to require fewer additional pointlike pion couplings. This 
tends to be true mostly for the Yang-Mills like versions. See for the analysis to order 
p 4 . In the sector involving Euvap this tends not to be so simple |89[ . 



5 Relation to other models 

As discussed in the introduction there are several variations on the theme of effective 
Lagrangians with quarks and mesons. In this section we describe how the ENJL model 



is related to the other approaches. This is an extended version of the discussion in |~3 



The relation with the Georgi-Manohar model and in particular the discussion about the 



pion-quark coupling can be found in f90|. 



For this comparison we first introduce a version that includes both bosonic and fermionic 
fields in the Lagrangian. Following the standard procedure of introducing auxiliary fields, 
we rearrange the Nambu-Jona-Lasinio cut-off version of the QCD Lagrangian in an equiva- 
lent Lagrangian which is only quadratic in the quark fields. For this purpose, we introduce 
three complex 3x3 auxiliary field matrices M(x), L^(x) and the so-called collective 

field variables, which under the chiral group G transform as 

M -> g R Mg{ (58) 

Lfi -> 9LL^g{ and -> g R R^g R . (59) 
We can then write the following identities: 

exp i J d 4 xCflj L (x) = 

r r ( N A 2 1 

J VM exp ij d 4 x I- (q L M^q R + h.c.) - —^-^tr(M^M)\ ; (60) 
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and 



expi J d 4 xC 
J VL^VR^ exp i J d 4 x 



NJL\ X ) 



N A 2 1 1 1 

™* L ^ + i^fej^J + < L - «)} • < 61 > 

where C^j L (x) and £^j L (:c) are the four-fermion Lagrangians in ( |T3D and QTSj) . 
By polar decomposition 

M = UH = £H£, (62) 

with {7 unitary and if (and H) Hermitian. From the transformation laws of M and £ in 
eqs. (0) and fl35|) , it follows that H transforms homogeneously, i.e. 

H ^h{%g L , R )Htf{Q,g LtR ). (63) 



The path integral measure in eq. (|60| ) can then also be written as 

exp i J d A xC s ^ L {x) = 

r r ( N A 2 | 

J ViVH exp ijd 4 xl- (q L £H£q R + q R ^q L ) - ^±jL-^H 2 \ . (64) 

We are interested in the effective action T e ff(H, £, L M , v, a, s,p) defined in terms of the 
new auxiliary fields H,^,L^, R^; and in the presence of the external field sources v^, a M ,s 
and p, i.e. 

e ir eff (H,t,L^;v,a,s, P ) = _ j eX p (-i J d^X-G^G^^j 

x | DgzDg exp i J d 4 x {qV QCD q + q L ^L^q L + q R ^R^ R - (q L ?H?q R + q R ^H^q L )} , 

(65) 

with Vqcd the QCD Dirac operator: 

Vqcd = Y{d„ + - *7 M K + 75fl M ) + *(* - H5P). (66) 

The integrand is now quadratic in the fermion fields. 

Here we can easily see how when we integrate out the quarks we will end up with 
different implementations of the vector fields. The fields and R^ correspond to the 
linear version discussed in the previous section. We can decouple the external fields and 

by doing a shift of the auxiliary vector fields 

L{R), - L'(i2% = L(R)^ + l{r), . (67) 
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Notice that this leads to precisely the type of mass term added in the gauged Yang-Mills 
vector description and the U and R' transform nonlinearly as gauge bosons under the 
chiral group. The relation with the CCWZ version will be given in section [8]. 

In principle we could also choose various versions for the scalars and pseudoscalars by 
the various choices possible for the matrix M. Two possibilities are shown in (|62|). H 
transforms in the CCWZ fashion [&J while H transforms as a purely lefthanded scalar. 

Most of the other quark-meson models described in the introduction are models con- 
taining quarks and pseudoscalars only. The QCD effective action model |J follows simply 
by setting 

= Rp = and H = (H) = M Q = -g s (qq) , (68) 

where (qq) is the quark vacuum expectation value derived in section |3|. The advantage of 
the present approach is that the spontaneous symmetry breaking that was added by hand 
in that model is now generated spontaneously. The approximations ( |B"8"D will be referred 
to later as the mean field approximations. 

The Georgi-Manohar model |6| requires a little more work to obtain. Here there is an 
additional free parameter, g^, the axial-coupling of the pseudoscalars to the constituent 
quarks. There have been some recent arguments about the order in N c this parameter is, 
see []90| and references therein. In the ENJL model it is obvious that this parameter is 
of leading order in 1/N C . In the purely fermionic picture it is obtained from the graphs 
shown in Fig. II. In general this parameter depends on the off-shellness of the pion but in 




Figure 4: The set of diagrams summed to obtain g^iQ 2 )- X is the insertion of the pion 
field and the other lines are fermions. 

the low-energy approximation it becomes a constant. In the language described above the 
parameter g^ appears due to the mixing of the pseudo-scalar fields and — R^. 

In general the quark-meson models include kinetic terms for the mesons as well. These 
are in the present approach of course assumed to be produced from the integration over 
the quarks. 
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6 Regularization independence 



The general method we will use to argue independence of the regularization procedure is 
the heat-kernel method. A review of this method can be found in . There exists various 
versions of the heat kernel method. The version we use here is the most naive one. More 



careful definitions also exist, see 45 and references therein. 



The underlying problem is that, as can be seen in Sect. |3], the chiral symmetry is 
spontaneously broken by the quadratic divergence. In a regulator that does not have the 
quadratic divergence, like dimensional regularization, one always works in the phase where 
chiral symmetry is explicitly realized in the spectrum. In the ENJL model this means 
that we treat it as being in a phase with weakly interacting massive quarks. The reason 
is that the logarithmic divergence in (|20|) has a negative sign so the vacuum energy from 
the logarithmic term is positive. To avoid this we have chosen a variation on the proper 
time regularization. Most regulators that preserve the presence of quadratic divergences 
do break the underlying chiral symmetry explicitly. The Ward identities have then to be 
used to determine the coefficients of the symmetry-breaking counterterms that have to be 
added to obtain chirally symmetric results. In general this is a very cumbersome method 
and we will use some simplified versions of it. 

In general we will consider several options. We can treat the heat kernel regularized 
by a specific regularization scheme. The one used here is the proper time heat kernel 
expansion. This is the scheme used to obtain the low-energy expansion of Sect. ||. We 
can then be more general in the heat-kernel expansion and leave the coefficients of the 
terms in the heat-kernel expansion completely free. This way we test a combination of the 
symmetry structure and the general couplings of the mesonic fields to the quarks only. It 
is rather surprising that in this case there are still several nontrivial results left. These 
type of results are in fact the major improvement of the methods used here as compared 



to the more traditional oneslJTO 



Since we would also like to go beyond the few first terms in the low-energy expansion 
it is necessary to either go to very high orders in the explicit heat kernel expansion or go 
to an alternative method where we directly regulate the Feynman diagrams. Here there 
are also several options. In it was shown how a regularization via dispersion relations 
and determining the subtraction constants from the heat kernel expansion can be used 
in this case. To go beyond two-point functions this method becomes very cumbersome 



as well and there a simpler method [|54|, |57|j was used. The essence of the method is to 
expand all one-loop diagrams of the constituent quarks into the basic integrals by removing 
all dependencies on the loop-momentum in the numerator via algebraic methods. All 
combinations that involve only Lorentz structures without g^ v are correctly reproduced this 
way. The Ward identities are then used to determine the Lorentz structures involving g^. 
For the two-point functions this procedure agrees with the dispersion relation technique and 
for 3 and higher point functions it agrees with the results from the heat-kernel expansion. 
The latter has been checked explicitly for the first few terms by comparing results from 
the full expansion with those from the heat kernel | f5"4"| . 

Let us now show the last procedure on the simplest example. We look at the one-loop 
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contribution to the two-point function 

IlJJ = i J d 4 xe tq - x (0\T (y* j (x)V kl (0)) |0> (69) 
with = q^j^qj- The relevant Feynman diagram is shown in Fig. |5|. We will here for 

O 



Figure 5: The one-loop fermion diagram. The lines are constituent quarks. The dots are 
insertions of the external currents. 

simplicity only quote the equal mass case. The resulting feynman integral expression after 
doing the Dirac algebra in four dimensions is proportional to 



/ 



d 4 x (-r 2 -q-r + m 2 )g^ + q^r u + r^q u + 2r^r u 
(27r) 4 (r 2 — m 2 )((q + r) 2 — m 2 ) 

This integral has to be proportional to 

QiiQu ~ q 2 9^ (71) 

from the Ward identities. Naively cutting of the integral in ( ffO|) leads to a piece of the 
form ([7T| ) but there is an extra term proportional to g^ v . This term should be absent and 
has to be removed via the Ward identities. 

We now remove the r 2 via r 2 = r 2 — m 2 + m 2 and q.r via 2q.r = ((q + r) 2 — m 2 ) — q 2 — 
(r 2 — m 2 ). This removes from the numerator a large fraction of the dependence on r. As 
the next step we combine the two numerators using a Feynman parameter 

i 1 
dx * . (72) 



(r 2 - m 2 )((q + r) 2 - m?) Jo (( r + xq) 2 - (to 2 - x(l - x)q 2 )) 2 

Then we perform in those integrals a shift to p = r + xq. The integral with an odd number 
of p's in the numerator vanishes. Those with two powers are of the form p^p v and are after 
integration proportional to g^ v . This procedure leads to the correct answer for the q^q v 
term but needs to have terms subtracted in the g^ v piece. This is done by requiring the 
full integral to be proportional to ([TT|). 
The final result is then proportional to 

/ o \ f 1 i f 74 X (l ~ X ) / \ 

[q^qu — q g^v) \ dx d p 2 . (73) 

v Jo J (p 2 — (m 2 — x(l — x)q 2 )) 
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The integral we now regularize via 



= / drr n - 1 exp (-tX) . (74) 

X n n-l)\Ji/M yy ' y ' 



1/AJ 

Rotating the p integral to Euclidean space finally leads to an answer proportional to 

m 2 — x(l — x)q 2 \ 




A 2 

x 



(75) 



with 



oo 



r(n,e) = / dTT n ~ Y e~ T . (76) 



This procedure can be easily generalized to the case with different masses and higher 
than two-point functions. The requirement of being proportional to (|7l|) is then replaced 
by using the appropriate Ward identities. 

The equivalent results to leaving the coefficients of the heat-kernel expansion free, is to 
find out which identities exist between the different one-loop Green functions and then to 
leave only the ones not related to others as completely free functions. In this case, similar 
to the low energy expansion, we are actually testing a whole class of models where the one- 
loop expressions are left completely free. A prominent example is the possible inclusion of 
extra low-energy gluonic effects as described earlier. 



7 The anomaly 

There have been claims, Fill and references therein, that the Extended Nambu-Jona- 



Lasinio model does not reproduce the correct QCD anomalous Ward identities. The correct 
result for the decay 7r° — ► 77 was found but there were deviations from the anomalous 
Ward identity prediction for the 77r°7r + 7r~ vertex. Here we review the solution of Ref. |58 



to this problem. A similar problem was encountered in constructing anomalous effective 
Lagrangians using full Vector Meson Dominance (VMD)P2"|j. The same solution also works 
in this case and it provides a simpler way to deal with the Ward identities than the 



subtraction method used in ref. p2 |. The point of view taken here is that the ENJL model 
is looked upon as a low energy approximation to QCD by only keeping the leading terms in 
1/A^. We know that the anomaly is a short-distance phenomenon that is not suppressed by 
the cut-off so these terms can be subtracted consistently to reproduce the correct anomalous 
Ward identities. The procedure here restores the correct terms. The lowest order terms 
thus become independent of the cut-off, but the higher order contributions (like 0(p 6 )) in 
the anomalous sector will still depend on the cut-off A x . 

We will first point out the underlying cause of the problem. This followed from the 
way the four quark vertices in |91] were treated. This is essentially equivalent to requir- 



ing VMD. The definition of the abnormal intrinsic parity part of the effective action for 
effective theories has already quite a history. After Fujikawa derived the anomalous Ward 
identities^] from the change in the measure in the functional integral Bardeen and 
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Zumino clarified the relation between the various forms of the anomaly found using this 
method p5| . This paper also clarified the relation between the covariant and the nonco- 
variant (or consistent) forms of the anomalous current. Leutwyler then showed how these 
different forms are visible in the definition of the determinant of the Dirac operator p6|. 



He also discussed the relation of the anomalous current to this determinant. At the same 
time Manohar and Moore showed how the Wess-Zumino term|43| can be derived from a 
change of variables in the functional integral in a constituent chiral quark model and how 
this can be used to relate different anomalously inequivalent effective theories [44]. 



What we will show here is that the terms that violate the anomaly generated by the 



procedure used in |9T[ can be subtracted consistently. We describe how the problem with 
the anomaly arises in the standard treatment of the ENJL model. Then we illustrate a 
simpler way to obtain the offending terms. This way will then show that these terms can be 
subtracted in a consistent fashion. We also show that our prescription does not influence 
the chirally covariant part of the effective action. 

Similar problems with the anomaly occur when one tries to formulate quark-meson 
effective Lagrangians which include vector and axial- vector meson couplings to the quarks. 
There the problem can be solved in a similar way by subtracting terms that contain only 
(axial)vector mesons and external fields. The same basic problem also occurs when trying 
to implement Vector Meson Dominance for the anomalous terms. We show how it is related 
to the problem in the ENJL model and can hence be solved similarly. Finally we explicitly 
state what our prescription corresponds to. 

In ( |i~2"|) the measure T>qD~q has to be defined in a way which reproduces the correct 



anomalous Ward identities. This means that the cut-off procedure should be defined with 
a Dirac operator that involves the external left and right handed vector fields. 

The standard way to analyze the generating functional flT2] ) is to introduce a set of 
auxiliary variables as described in section [| to obtain an action bilinear in fermion fields. 
We will concentrate here on the vector-axial- vector part since it is that one that may 
generate the problems with the anomaly. The scalar-pseudoscalar part is already treated 
in ref. |pf[ . 

Formally the Lagrangian in the exponential can be rewritten in terms of the full fields 
I'fi = ^ + r 'fi = r v + ^ an d s \v' ■ The latter are defined by s' — ip'^s = s — ip^ 5 + 
M"f L + M^ R . We then have that 

r(l,r,s,p,L,R,M) = T(l'y,s',p') . (77) 

with r(/,r, s,p) defined by 

expir(Z, r, s,p) = J VqVqexp J d A x (q^iD^q — q (s — ipjs) g)| . (78) 

We can then integrate out the fermions to obtain the effective generating functional 
as a function of the external fields and the auxiliary fields. There is one caveat here and 
that is precisely the cause of the problem observed in |91||. The measure that corresponds 



to the standard procedure is then defined by a Dirac operator that is a function of l'ni r 'u, 
rather than a function of and r M . 
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Let us show in a compact fashion how this problem occurs. For simplicity we temporar- 
ily neglect the scalar-pseudoscalar part. The effective action T(l',r',s,p) can be related 
simply to T(l, r, s, p) by introducing the fields 



and 



(79) 



These fields transform under the chiral symmetry group in the same way as l^, r^. We can 
now describe the effective action after integrating the fermions as 

r(vy,s, P ) = r(i,r,s,p) + £dt-r(i t ,r t ,s lP ) 

= T(l,r,s,p) + £dttv(L^r(l\r\s,p)+R^r(l t ,r t ,s,p)\ .(80) 



The last two terms in (|80| ) correspond to the left and right handed current. This current 
consists out of two pieces, a non-anomalous and an anomalous part. The part that is non- 
anomalous causes no problem and one can use the standard heat kernel methods as used 



in Refs. |T0| , 13], section to obtain information about the generating functional fll2|) . 

The anomalous part of the current can also be written as the sum of a local chirally 
covariant part and a local polynomial of 0(p 3 ) in /*, r*[95]. If we now insist that at the 
first step, where we integrate out the fermions, we should have the global chiral symmetry 
exact (this corresponds to choosing the left-right form of the anomalous current) this local 
polynomial contains two pieces. One is a function of Z* and its derivatives only and the 
other one is a function of r* and its derivatives. This globally invariant form is precisely 
the form that a "naive" application of the heat kernel method would give ||96|| . 

The anomalous left and right currents [of 0(p 3 )] in eq. fl30| ) have the following form in 
the left-right symmetric scheme 



6T(l,r,s,p) 



SI" 

ST(l,r,s,p) 



jR I -ti . 

J H * J ii i 



R . 



N, 



487T 2 

-^^[il v l a l p + {l vw lp}\ ; 

d ll l{r) v - dj(r)n - i[l(r)^ l(r) u ] 
iU^dM-ir^U + iUL) . 



The matrix U is the "phase" of M. U = £ 2 with M = C.HC,. Here H is hermitian 
and £ is unitary. The currents and jj? can be obtained from and jjj by a parity 
transformation. 

Since in (|80D the part T(l,r,s,p) already saturates the inhomogeneous part of the 
anomalous Ward identities the remainder should be locally chirally invariant. The parts 
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that are not locally invariant in the last two terms of eq. (|80| ) should thus be subtracted. 
As can be seen from (|8l|) these terms are a local function of /, L and their derivatives (plus 
the right handed counterpart). The change in the definition of the measure involves only 
the fields Z, L, r and R so the local terms that can be added to the effective action to obtain 
the correct Ward identities should only be functions of these and their derivatives. The 
preceding discussion shows that the terms that spoil the anomalous Ward identities are 
precisely of this type. 

As a consistency check we will show that the contribution of the local chirally covariant 
part of the anomalous current to the resulting effective action can not be changed by adding 
globally invariant counterterms that are functions of I, L and their derivatives only. The 
full list of terms that could contribute is (an overall factor of e^ VOL ^ is understood). 

^{L^L u L a Lp) , tr(D fl L u L Q L/ 3 ) , 

tr(l^L a Lp) , tiil^DaLp) (82) 

and their right handed counterparts. All others are related to these via partial integrations. 
The first term vanishes because of the cyclicity of the trace. The second one is a total 
derivative. The third one is forbidden by CP invariance and the last one vanishes because 
of the Bianchi identities for l^ u . 

This is just proving that the standard procedure of adding counterterms and deter- 
mining their finite parts by making the final effective action satisfy the (anomalous) Ward 
identities also works here giving an unambiguous answer. 

We have used the left-right symmetric form of the anomaly. But it is obvious from 
the discussion above, that by following an analogous procedure to the one given here 
in any scheme of regularization of the chiral anomaly one obtains the same result since 
the difference between the anomalous current in two of these schemes is a set of local 



polynomials that can only depend on l'^ and r^, ref. [p6 |. A scheme of particular interest 
is that where the vector symmetry is explicitly conserved. In order to obtain this form of 
the Wess-Zumino action one has to add a set of local polynomials that only depend on V 
and r' to the left-right symmetric one. These are given explicitly in ref. |)3| . 



In the basis of fields we have been working until now the effective action in the non- 
anomalous sector has generated a quadratic form mixing the pseudoscalar field and the 
axial- vector auxiliary field (roughly speaking — L^). It is of common practice to change 
to a basis where this quadratic form is diagonal (e.g. see |13[]). Afterwards the vector 
and axial-vector degrees of freedom can be removed by using their equations of motion to 
obtain an effective action for the pseudoscalars only. In this way one also introduces the 
axial coupling, the so-called gA, which in the chiral constituent quark model || corresponds 
to the axial vector coupling of constituent quarks to pseudoscalar mesons. In our effective 
action, this change of basis can only generate local chiral invariant terms that therefore 
cannot modify the Wess-Zumino effective action [|^] and the standard predictions at order 
p 4 for 7T° — > 77 and 77r + 7r~7r° will be satisfied. There will of course be changes at higher 
orders due to the chiral local invariant terms. Thus, the value of gA is not constrained by 
the chiral anomaly which is a low-energy theorem of QCD contrary to the conclusion of 
ref. |JT| and in agreement with the results of ref. |H] . 
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These changes due to higher orders are very similar to the description using the hidden 
symmetry approach |)7| (see also [p8fl ) and the gauged Yang-Mills approach as given in 
ref. [39]. This prescription is also precisely the prescription that was used in ref. [61] 
to construct the lowest order anomalous effective chiral Lagrangian involving vector and 
axial-vector fields and obtain predictions for the "anomalous" decays of these particles 
within the ENJL model. 

|T 



We would like to add one small remark about ref. 



In this reference only the p 



equations of motion for L^, were used. In principle there is also a contribution from 
the p 4 part proportional to e^ VOi " when substituted into the p 2 term of the effective action. 
This contribution does however cancel between the p 2 term and the "mass" term for the 
auxiliary fields. 

Our discussion was in the framework of the ENJL model. The root of the problem was 
the relation ([T7]). As mentioned before similar problems occur in effective quark-meson 
models with explicit spin-1 mesons couplings to the quarks and in the old approaches that 
require full vector meson dominance (VMD). The basic requirement of VMD is that vector 
mesons couple like the external fields. If we describe the physical vector mesons by fields 
L M and R^, this requirement can be cast in the form {(f) stands for all the other fields 
involved) 



5 n T(l,r,L,R, 



d'-'L^S" 



R v 



5 n T(l,r,L,R, 



L=R=0 



L=R=0 



for < m < n . 

Using the Taylor expansion of T(7, r, L, R, 0) in /, r, L and R and applying eq. 
be shown that the action then only depends on I' = I + L and r' = r + R, i.e. 



T(l,r,L,R, 



T(l + L,r + #,0,0, 



(83) 
it can 

(84) 



This will lead to precisely the same type of problems as seen in the ENJL model since this 
is the same relation as eq. flT7|). Here again it can be remedied by adding local polynomials 
in l,L,r and R precisely as was done before. 

Now, what does our prescription mean in practice? It means that vector and axial- 
vector fields are consistently introduced in the low-energy effective Lagrangian by requiring 
a slightly modified VMD relation 



S n T(l, r, L, R, > 



Rv 



5 n T(l,r,L,R, 



L=R=0 

for < m < n with m + n > 1 



S m l ll S n ~ m r 



L=R=0/ 



local 
covariant 



J5) 



instead of the usual VMD requirement in eq. (|83|). This is equivalent to use the standard 
heat kernel expansion technique (for a review see |45|) for the non-anomalous part, i.e. no 
Levi-Civita symbol, and for the chiral orders larger than p 4 in the anomalous part, i.e. 
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terms with a Levi-Civita symbol. For the 0(p 4 ) part of the anomalous action one has just 
the usual Wess-Zumino term and for the chiral orders smaller than p 4 in the anomalous 
part one has 

J dt tifr tf(l\r\U) + R*J*(l\r\U)) . (86) 

Here the anomalous currents J^' R are those defined in eq. (81). 

In the present work we have been implicitly using a representation similar to the so- 
called vector model (model II in ref. [|55|] ) to represent vector and axial- vector fields as the 
most natural way within the ENJL model we are working with. However, it is straightfor- 
ward to work out the analogous prescription to eq. (B3T) for any other suitable representa- 



tion of vector and axial- vector fields (tensor, gauge fields, ■ • -)[|85j to implement VMD in 



both the anomalous and the non-anomalous sectors of the effective action. 



8 Analysis to order p A 

8.1 The Mean Field Approximation 

We shall first discuss a particular case of T e ff(H, £, L M , R^; v, a, s,p) as defined in eq. (|65|) . 
It is the case corresponding to the mean field approximation, where 



H{x) =< H >= MqI, 



and 



R„ = 0. 



J7) 



The effective action T e ff(MQ, £, 0, 0; v , a, s,p) coincides then with the one calculated in 
Ref. ||, except that the regularization of the UV behaviour is different. In Ref. |§, the 
regularization which is used is the ( function regularization. The results, to a first approx- 
imation where low-frequency gluonic terms are ignored, are as follows: 



Nr. . , r2r - Mr 



167T 



;4M^r(0, 



A 2 

x 



(88) 



and 



fo B o 



< tptp >- 



16tt 2 q v ' A 2 ' 



for the lowest 0(p 2 ) couplings of the low-energy effective Lagrangian in 
For the 0(p 4 ) couplings which exist in the chiral limit we find 



39) 



2L X 



N c 1 w M% 

-r(2, 



16vr 2 12 



A 2 

x 



N c 1 
16tt 2 6 



rri. 



A x 



- 2r(2, 



A x' 



(90) 
(91) 
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for the four derivative terms; and 



N r 1. 



M, 



Q 



16tt 2 3 



N r 1 



'10 



167T 2 6 



rHl, 



(92) 



(93) 



for the two couplings involving external fields. If one lets Mq/A 2 — > 0, then T(n, 0) = 
r(n) = (n — 1)! for n > 1, and these results coincide with those previously obtained in 
Refs. |7|, H, ||, § and |4§ to |9|. 



When terms proportional to the quark mass matrix Ai are kept, there appear four new 
Lj couplings (see eq. (|32|)). With 



Mr 



(Sol | < ipij) > 
the results we find for these new couplings are 



N c p 
16tt 2 2 



r(o, 



L 7 



N c 1 
" 16tt 2 12 

Nr. 1 



-pT(0, 



= 

= 

A x 



rn, 



16vr 2 24 



+ i r ( i, 

Mr 



A 2 

x 



6p( P -i)r(o,-^) + r(i 



A l 



(94) 

(95) 

(96) 
(97) 
(98) 

(99) 



If we identify T(0, Mq/A 2 ) = \og(p 2 /M^), and take the limit T(n > 1, il^/A 2 -»■ 0), these 
results coincide then with those obtained in Ref. ||. (Notice that p is twice the parameter 
x of Ref. p.) 

The fact that L4 = L$ = and L2 = 2Li is more general than the model calculations 



we are discussing. As first noticed by Gasser and Leutwyler |69| , these are properties of the 
large- N c limit. The contribution we find for L 7 is in fact non-leading in the 1/N C expansion. 
The above result is entirely due to the use of the lowest-order equations of motion (see 
the erratum to Ref. ||). In the presence of the £/a(1) anomaly, L 7 picks up a contribution 
from the rf pole and becomes 0(N 2 ), |69| []. 

Finally, we shall also give the results for the Hi and H 2 coupling constants of terms 
which only involve external fields: 



N c 1 
16tt 2 12 



2r(0, 



Ml 
A 1 



-rri. 



Ml 
A x 



(100) 



4 This counting is somewhat misleading since it first relies on l/N c to be small to have ml, of order 
1/N C and then expands in 1/m^,, or 1/N C large. 
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16vr 2 12 



1 /- /bf2 W-' 1 

6p 2 r(-i,^)-6p(p + i)r(o,^) + r(i,^ 

X X X J 



(101) 



8.2 Beyond the Mean Field Approximation 

In full generality, 

H(x) = M Q 1 + a(x) , (102) 

and the effective action T e ff(H,£, L M , R^\v,a, s,p) has a non-trivial dependence on the 
auxiliary field variables cr(x), L^x) and R^{x). It is convenient to trade the auxiliary left 
and right vector field variables L M and R^, which were introduced in eq. (|6lD , by the new 
vector fields 

=£L^±^R^. (103) 

From the transformation properties in eqs. (|35|) and fl5U|), it follows that transform 
homogeneously, i.e. 

WM^h($,g)W^tf($,g). (104) 

We also find it convenient to rewrite the effective action in eq. ( p5| ) in a basis of constituent 
chiral quark fields 

Q = Ql + Qr and Q = Q L + Q R , (105) 

where 

Ql = iqh , Ql = qL$ ] Qr = ^Qr , Qr = QrL (106) 
which under the chiral group G, transform like 

Q -> h($, g)Q and Q -> Q/i($, (107) 

In this basis, the linear terms (in the auxiliary field variables) in the r.h.s. of eq. ( |65] ) 
become 

Q (-H + i 7 ^< +) - ^ _) ) (108) 

At this stage, it is worth pointing out a formal symmetry which is useful to check 
explicit calculations. We can redefine the external vector-field sources via 

Jm-?, = Jm + ^ ( 109 ) 

r^r'^r^ + R, (110) 

and 

M^M'{x) = M + £a{x% (Ill) 

The Dirac operator T>e, when reexpressed in terms of the "primed" external fields, is for- 
mally the same Dirac operator as the one corresponding to the "mean field approximation." 
In practice, it means that once we have evaluated the formal effective action 

exp r ef f(Ap,M) = J VQVQexp J d A xQV E Q = detV E (112) 

we can easily get the new terms involving the new auxiliary fields L M , i? M and a by doing the 
appropriate shifts. The formal evaluation of T e ff(A^, M) to 0(p 4 ) in the chiral expansion 
has been made by several authors (see Refs. [|4TJ] to 
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8.3 The constant g& and resonance masses 

When computing the effective action V e ^{A^ M) in eq. ( |112| ), there appears a mixing 
term proportional to Xx^W^ . More precisely, one finds a quadratic form in £^ and 
(in Minkowski space-time): 



with 



and 



a < W { ~\W { ~^ >+(3< ^W { -^ > +7 < > (H3) 
N r /1A? Mr 



a 



16tt 2 



AT Af 1 

/5 = - 2 T^ r (0, ^)MJ and 7 = . (115) 

W H „ -> WM + (1 - , (116) 

diagonalizes the quadratic form. There is a very interesting physical effect due to this 
diagonalization, which is that it redefines the coupling of the constituent chiral quarks to 
the pseudoscalars. Indeed, the covariant derivative in eq. (|112|) becomes 



The field redefinition 
with 



V E = 7 ^V M = 7 „ (fy + iGp + T, - ~^(gA^ - WV) - V«J . (118) 

qa can be identified with the qa coupling constant of the constituent chiral quark model 
of Manohar and Georgi || . 

In the calculation of T e fj(A^, M) we also encounter kinetic-like terms for the fields 
and W^ + \ Comparison with the standard vector and axial- vector kinetic terms requires a 
scale redefinition of the fields and to obtain the correct kinetic couplings, i.e. 

V, = \ V W^, At = \ A Wjr\ (119) 

with 



and 

16vr 2 3 



, M o, , M o 

r(o,-rf) - r(i, -J- 

X X 



'12V, 



These and fields are the ones that transform in the standard CCWZ way [Q for 
the vector and axial-vector fields. 
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This scale redefinition gives rise to mass terms (in Minkowski space-time) 



l ~M 2 v tT{V»V») + -Mitr(^) 



with 



My 



and Ml 



X 2 



2a 



'1221 



(123) 



The same comparison between the calculated kinetic and mass terms in the scalar 
sector, with the standard scalar Lagrangian in eq. (|4"9"D, requires the scale redefinition 



with 



\2 



N c 2 
16tt 2 3 



S(x) = \ s a(x), 
M 2 

3r(o,^#) -2rd, 



A 2 

x 



The scalar mass is then 



N c 8M 2 / M 2 
Ms ~ 167T 2 A 2 (U ' 



s 



A 2 

X 



(124) 
(125) 
(126) 



8.4 The couplings of the C^jj Lagrangian 



The Lagrangian in question is the one that we have written in section [|, in eqs. (49), 
( |50D and (|5T|), based on chiral-symmetry requirements alone. These requirements did not 
fix, however, the masses and the interaction couplings with the pseudoscalar fields and 
external fields. The results for the masses which we now find in the extended Nambu- 
Jona-Lasinio model are given by eqs. fll23| ) and (|126|) in the previous subsection. These 
are the results in the limit where low-frequency gluonic interactions in cIqcd i n ec L- (0) are 
neglected, i.e. the results corresponding to the first alternative scenario we discussed in the 
introduction. For the other coupling constants, and also in the limit where low-frequency 
gluonic interactions are neglected, the results are: 



Mr 



4*=16^ r <°'Af>' 
instead of the mean field approximation result in eq. (^8|):p| 

fv = V2X V , f A = V2g A X A 

N c 1 y/2 



9v 



16tt 2 Av 6 



;i-9l)r(o ; 



2^ir(i, 



M, 



2 1 



A 2 

X J 



(127) 

(128) 
(129) 



5 This implicitly changes the value of B via eq. (R9 
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for the vector and axial- vector coupling constants in (|50| ) and ([51]); and 



N c M Q 

167T 2 Ac 



p 



rr-i. 



A x 



2r(o, 



A x 



r(o, 



Ax 



for the scalar coupling constants in (|49|) . 

There are a series of interesting relations between these results: 



Mi 



3 A 2 
2 _ 6 IV X 



v 



(130) 



131} 



(132) 



with the two solutions 



and 



>, r(i,^) 

M A 1 



2a 



+ 6M 2 Q , 



fyMy 2 



^ = and = 



fv^v = flMl + fl 



(133) 

(134) 

(135) 
(136) 



The last relation is the first Weinberg sum rule [ 100| . Using this sum rule and the second 
solution for gA, we also have 

9a — 1 -J^- (137) 



mi 



Therefore gA < 1. The two relations in eqs. ( |136j ) and ( |137| ) remain valid in the presence 
of gluonic interactions, i.e. the gluonic corrections do modify the explicit form of the 
calculation we have made of fv, My and gA, but they do it in such a way that eqs. ( |136| ) 



and (|137|) remain unchanged. 



8.5 The coupling constants Lj's, Hi and H2 beyond the mean field 
approximation 

These coupling constants are now modified because we no longer have gA — 1- With the 
short-hand notation 



x 



(138) 
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the analytic expressions we find from the quark-loop integration are the following: 



= 2Li = 

N c 1 
16tt 2 24 



N c 1 
16tt 2 24 
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^4 = 0, 

AL n , 
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16tt 2 12 
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6p 2 r(-l, x) - 6p(p + ^)r(0, x) + g A F{l, x) 



(139) 

(140) 

(141) 

(142) 

(143) 
(144) 

(145) 
(146) 
(147) 
(148) 
(149) 



3, 5 and 8) as well as Hi receive explicit contributions 



from the integration of scalar fields. This is why we write 



3, 5, 8; 



H2 = H2 + H2 with Li, H2 the contribution from the quark-loop and Lf , H2 that from 
the scalar field. The results for L\, L2 and L3 agree with those of Ref. [ 101 ], where these 
couplings were obtained by integrating out the constituent quark fields in the model of 
Manohar and Georgi ||. At the level where possible gluonic corrections are neglected, 
the two calculations are formally equivalent. The results for L 4 to L w agree with those of 



Ref. [37 



We note that between these results for the Li's, Hi and the results for couplings and 
masses of the vector and axial- vector Lagrangians, which we obtained before, there are the 
following interesting relations: 



no 



f 2 A ) and 2H X 



-\(fv + f. 



(150) 



(151) 



As we shall see in the next subsection, these relations, like those in eqs. ( |136| ) and ( |137| ), 
are also valid in the presence of gluonic interactions. The alerted reader will recognize that 
these relations are precisely the QCD short- distance constraints which, as discussed in 
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Ref. |85| , are required to remove the ambiguities in the context of chiral perturbation 
theory to 0(p 4 ) when vector and axial- vector degrees of freedom are integrated out. They 
are the relations which follow from demanding consistency between the low-energy effective 
action of vector and axial-vector mesons and the QCD short-distance behaviour of two- 
point and three-point functions. It is rather remarkable that the simple ENJL model we 
have been discussing incorporates these constraints automatically. 

There is a further constraint that was also invoked in Ref. f85|. It has to do with the 
asymptotic behaviour of the elastic meson-meson scattering, which in QCD is expected 
to satisfy the Froissart bound ||102|| . If that is the case, the authors of Ref. |85j concluded 
that, besides the constraints already discussed, one also must have 



;9v 



2Li ; L 3 



-6L a . 



(152) 



As already mentioned, the second constraint is a property of QCD in the large N c limit. 
The first and third constraints, however, are highly non-trivial. We observe that, to the 
extent that 0(N c g A ) terms can be neglected, these constraints are then also satisfied in 
the ENJL model. 



When the massive scalar field is integrated out |84j, there is a further contribution to 
the constants L3, L5, L§ and H2 with the results: 
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(153) 

(154) 

(155) 
(156) 



This result for Lf disagrees with the one found in Ref. f37|. Also, contrary to what is 



found in Ref. pYfl , there is no contribution from scalar exchange to L 2 . 

It is interesting to point out that L 5 , Lf and L 8 , Lf each depend explicitly on the 
parameter p. This dependence, however, disappears in the sums 



Nc 1 3 



L5 + L *- 16VM 



and 



L 8 + La 



1 f 2 9 A 



g\ [r(0,x)-r(l,a;)] 
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(157) 
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Similar simplified expressions for the other Lj can be found in |64] 
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8.6 Results in the presence of gluonic interactions 

The purpose of this section is to explore in more detail the second alternative, which we 
described in the introduction, whereby the four-quark operator terms in eqs. fllif) and fllBD 
are viewed as the leading result of a first-step renormalization a la Wilson, once the quark 
and gluon degrees of freedom have been integrated out down to a scale A x . Within this 
alternative, one is still left with a fermionic determinant, which has to be evaluated in 
the presence of gluonic interactions due to fluctuations below the A x scale. The net effect 
of these long-distance gluonic interactions is to modify the various incomplete gamma 
functions T(n,x = Mq/A 2 ), which modulate the calculation of the fermionic determinant 
in the previous sections, into new (a priori incalculable) constants. We examine first how 
many independent unknown constants can appear at most. Then, following the approach 
developed in Ref. |§, we shall proceed to an approximate calculation of the new constants 
to order asN c . 



8.6.1 Book-keeping of (a priori) unknown constants 

The calculation of the effective action in the previous sections was organized as a power 
series in proper time. 

In the presence of a gluonic background, each term in the effective action, which orig- 
inates on a fixed power of the proper-time expansion of the heat kernel, now becomes 
modulated by an infinite series in powers of colour-singlet gauge-invariant combinations 
of gluon field operators. Eventually, we have to take the statistical gluonic average over 
each of these series. In practice, each different average becomes an unknown constant. If 
we limit ourselves to terms in the effective action to 0{p i ) at most, there can only appear 
a finite number of these unknown constants. We can make their book-keeping by tracing 
back all the possible different types of terms that can appear. 

In the presence of gluonic interactions, there then appear 10 unknown constants: 7_i; 
7oi, 702, 703; 7ii) Ti2, Ti3, 7i45 7n, 722- To these, we have to add the original G s and G v 
constants, as well as the scale A x . However, the unknown constant (1 + 7_i) in eq. ( |161[ ) 
can be traded by an appropriate change of the scale A x , 

M 2 

r(-l,x)=r(-l,z){l+ 7 _i}; x = ^f, (159) 



A 2 

x 



and a renormalization of the constant Gs, 

A 



Gs-^G S = -f G s . (160) 

x 



Altogether, we then have 12 (a priori unknown) theoretical constants and one scale A x . 
They determine 18 non-trivial physical couplings (in the large- N c limit) of the low-energy 
QCD effective Lagrangian: < Tpip >, f w , L u L 3 , L 5 , L 8 , L 9 , L 10 , H u H 2 , f v , /a, 9v, c, 
c d , M s , M v and M A . 



Ill ; 
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In full generality, the results are: 



< ^ >= 



;4M3r_ 1 (l + 7 -i), 



16tt 2 q 

167r2 M ^r (i + 7 oi); 



I f 2 = N c S/f 2 



(161) 
(162) 



L 2 — 2L\ 

,2/i 2- 



N c 1 

16tt 2 24 X 
3 I 



(163) 



(i - <?A) 2 r (i + 7o3 ) + %i(i - sDiyi + -7i2 - -713) + 2^r 2 (i + 721 ) 



N c 1 
16tt 2 24 



-3(1 - ^i) 2 r (i + 7o 3 ) + 4^r x (i + 713) 



-i2^(i - ^)ri(i + ^712 - - 8<?lr 2 (i + 1( 721 + 722 )) 



N c 1 3 I+701 



' " 2M 2 ' 



L5 ^4^rr^^ 1+ ^- r ^ 1+ ^' 



Lo = 



N c 

" 167T 2 
N n 1 



1 1 + 701 1 Ti(l +713) 



#A r o(l + 7oi), 



16tt 2 6 
-^10 — — 



16 1+702 24r (l+7oi) 
(1 - g 2 A )T (l + 703 ) + 2g 2 A T 1 (l + ^ 7 i 2 - ^ 7l3 ) 

A< ' (l-^ 2 )r (l + 703) +^(1+713) 



#1 



167T 2 6 

N c 1 
16tt 2 12 



(i + ^)r (i + 7o3 ) -5^1(1 + 713) 



(164) 
(165) 
(166) 
(167) 

(168) 
(169) 
(170) 



Ho 



N c 1 
16tt 2 12 



6p 2 T^(l + 7 _0 - 6p 2 r (l + 7o2 ) - 6p^r (l + 701) + tfUMl + 713) 



and 
with 



rri m . 

2 " M 2 ' 
/v = V2\ v 

f A = V2g A X A , 

^ = !&^ r ° (1 + 703) 
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(171) 
(172) 

(173) 

(174) 

(175) 



and 



with 



9v 



^A = ^- 2 \W+103)-m+ ll3 )}. 
l07T z 6 

i^ x^ir [ (1 ~ ^ )r ° (1 + 7os) + 2 ^ ri(1 + ^ 7l2 " h n) 

A " % [r^i + 7 _ 1 )-2r (i + 702 )], 



Cd 



167T 2 A 5 

M Qn 2 



167T 2 A< 



2^1 [r„(l + 7oi)-ri(l + T ii) 



N 2 

■ " -[3r (l + 7oi)-2r 1 (l+7 14 )]; 



(176) 
(177) 
(178) 
(179) 

(180) 



M 



2 _3 A* 



2G y (A 2 )r(0,x)(l + 7 03) 



Mi 



1 - 



r(i,a0(i + 7i 3 ) 



M^ + 6Mji±^. 

1 + 703 



181) 
182) 
183) 



r(0,a;)(l + 7o3), 

There exist relations among the above physical couplings which are independent of the 
unknown gluonic constants. They are clean tests of the basic assumption that the low- 
energy effective action of QCD follows from an ENJL Lagrangian of the type considered 
here. The relations are 



/£M£ - fZM% = f; (first Weinberg sum rule), 



184) 



U = \fvgv , (185) 
Lxo = + \f A , (186) 

2#! = -\fy - \fl , (187) 

and 

H ? + 2L 8 c m 

= ~c~d ■ (188) 

The first four relations have already been discussed in the previous subsection. The 
combination of couplings in the r.h.s. of eq. (|188|) is the one that appears in the context of 
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non-leptonic weak interactions, when one considers weak decays such as K — > ttH (light 
Higgs) ||103 |. In fact, from the low-energy theorem derived in |69| it follows that 



1 <0|uu|0> 



<0|ss|0> _ ^ 

(189) 



2L 5 4 f K /U - 1 ' 
Experimentally 

f K /U~l = 0.22 ± 0.01 . (190) 

Unfortunately, the numerator in the r.h.s. of (|189|) is poorly known. If we vary the 
ratio 

<ss> _ 1 f rom _o.l to -0.2, (191) 

< uu > 

as suggested by the authors of ref. ||103|| , then eq. (|188|) leads to the estimate 



Cm / Cd = -1.1 x 10" 1 to -2.3X10" 1 . (192) 

With this estimate incorporated in eq. (|56|), we are led to the conclusion that 

\c d \ ~ 34 MeV . (193) 

In the version corresponding to the first alternative, the results for c m and q are those 
in eqs. ( |130| ) and (|131|) . We observe that in this case c m /cd comes out always positive for 
reasonable values of Mq/A"^. In fact from the gap-equation discussed in section |3| it is 
obvious that (qq) increases with increasing current quark mass for not too high current 
masses. 

8.6.2 Gluonic correction to 0(asN c ) 

We can make an estimate of the ten constants 7_i; 701, 702, 103', 7n, 712, 7i3, Ju] 721 and 
722, by keeping only the leading contribution, which involves the gluon vacuum condensate 
< ^-GG > /Mq as was done in Ref. ||. The relevant dimensionless parameter is 

2 < OSQQ > 

(194) 



^ 6N C Mq 

Notice that in the large-iV c limit, g is a parameter of 0(1). One should also keep in 
mind that the gluon average in (|194|) is the one corresponding to fluctuations below the 
A x scale. The relation of g to the conventional gluon condensate that appears in the QCD 
sum rules [|J |104|| is rather unclear. We are forced to consider g as a free parameter. Up 
to order 0(asN c ), this is the only unknown quantity which appears, and we can express 
all the 7's in terms of g. We find: 

T(l,x) , \ 

^ = rpti 29; < 195 > 
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r(4,x) / o 



(196) 



(197) 



(198) 



" r(2,x) V2/5, 

Notice that the combination ~7i 2 — 5713 entering some of the Lj's coupling constants 
is zero. This is the reason why it was found, in Ref. ||, that in the limit — > 1, L 2 and 
L 9 have no gluon correction of 0(a s N c ). 

To this approximation, we have then reduced the theoretical parameters to three un- 
known constants Gs, Gy and g, and the scale A x . 



8.7 Discussion of numerical results 

In the ENJL model, we have three input parameters: 

G s , G v and A x . (199) 
The gap equation introduces a constituent chiral quark mass parameter Mq, and the ratio 



x 



is constrained to satisfy the equation 

7 ^- = xr(-i,x)(i + 7 _ 1 ). 

Once x is fixed, the constants g& and Gy are related by the equation 

1 

9 A 



l + 4G y xr(0,x) (I+701) 



(200) 



(201) 



(202) 



Therefore, we can trade Gs and Gy by x and g&\ but we need an observable to fix the 
scale A x . This is the scale which determines the p mass in eq. (|182|) , i.e. 



2 M 2 
Al = -M v G v m^-)(l + l03 ). 

x 



(203) 



There are various ways one can proceed. We find it useful to fix as input variables the 
values of Mq, A x and g&. Then we have predictions for 



(204) 
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M s , M v and M A (205) 

fv , g v and f A ; and c m and c d (206) 

and the 0(p 4 ) couplings: 

Li (i = 1, 2, . . . , 10) and #i , ff 2 . (207) 

In principle we can also calculate any higher-Q(p 6 ) [ |105|1 coupling which may become of 
interest. So far, we have fixed twenty- two parameters. Eighteen of them are experimentally 
known. 

In the first column of Table [I] we have listed the experimental values of the parameters 
which we consider. In comparing with the predictions of the ENJL model, it should be 
kept in mind that the relations (|184|) to (|186|). are satisfied by the model while 



Li = \&\ L 2 = 2L 1 = ^g v ; L 3 = -3L 2 (208) 

only have numerically small corrections. These relations are rather well satisfied by the 
experimental values and thus constitute a large part of the numerical success of the model. 

We have also used the predictions leading in 1/N C , so that we have L\ = L 2 /2, L4 = 
Lq = 0, and we do not consider L7 since this is given mainly by the if contribution ||69| . 
In evaluating the predictions given in Table ^ we have used the full expressions for the 
incomplete gamma functions and the numerical value of the 7^ in terms of g given in 
eqs. QB^) to (PSp. 

The first column of errors in Table |I| shows the experimental ones. The second column 
gives the errors we have used for the fits. When no error is indicated in this column, 
it means that we never use the corresponding parameter for fitting. This is the case for 
< qq >, which is quadratically divergent in the cut-off and which is not very well known 
experimentally. This is also the case for c m , which depends on < qq >. Fit 1 corresponds 
to a least-squares fit with the maximal set of parameters and requiring g > 0. Fit 2 
corresponds to a fit where only /„- and the Li are used as input in the fit, while fit 3 has 
the vector and scalar mass as additional input. The next column, fit 4, is the one where 
we require g A — 1, i-e. we start with a model without the vector four-quark interaction. 
Here there are no explicit vector (axial) degrees of freedom, so those have been dropped 
in this case. This fit includes all parameters except My, M A , fv, gv and f A . Finally, fit 
5 is the fit to all data, keeping the gluonic parameter g fixed at a value of 0.5. The main 
difference with fit 1 is a decrease in the value of Mq. The value of A x changes very little. 
In addition the result with the constraint Gs = 4Gy, ([Rf), included is shown as fit 6. Fit 
7 is the result without gluonic corrections and Gy = as suggested by JFD . 



The expected value for the parameter g, if we take typical values from, e.g, QCD sum 
rules, is of 0(1). None of the fits here really makes a qualitative difference between a g of 
about 0.5 to 0. Numerically we can thus not decide between the two alternatives mentioned 
in the introduction. This can be easily seen by comparing fit 1 and fit 5, or fit 4 and fit 7, 
in Table |. 
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In all cases acceptable predictions for all relevant parameters are possible. The scalar 
sector parameters tend all to be a bit on the low side; but so is the constituent quark 
mass. The predictions for the Lj's are reasonably stable versus a variation of the input 
parameters. For L5 and Lg, this is a major improvement as compared with the predictions 
of the mean field approximation || . The typical variation with input parameters can be 
seen in table 2 of [O . 





exp. 


exp. 


fit 


fit 1 


fit 2 


fit 3 


fit 4 


fit 5 


fit 6 


fit 7 




value 


error 


error 
















f 

Jit 






10 


89 


86 


86 


87 


83 


86 


86 


v ^ yy ^ 








981 


960 


255 


1 78 




210 


1 70 


10 3 • L 2 


1.2 


0.4 


0.5 


1.7 


1.6 


1.6 


1.6 


1.7 


1.5 


1.6 


10 3 ■ L 3 


-3.6 


1.3 


1.3 


-4.2 


-4.1 


-4.4 


-5.3 


-4.7 


-3.1 


-3.0 


10 3 • L 5 


1.4 


0.5 


0.5 


1.6 


1.5 


1.1 


1.7 


1.6 


2.1 


1.9 


10 3 • L 8 


0.9 


0.3 


0.5 


0.8 


0.8 


0.7 


1.1 


1.0 


0.9 


0.8 


10 3 • Lg 


6.9 


0.7 


0.7 


7.1 


6.7 


6.6 


5.8 


7.1 


5.7 


5.2 


10 3 ■ L10 


-5.5 


0.7 


0.7 


-5.9 


-5.5 


-5.8 


-5.1 


-6.6 


-3.9 


-2.6 


10 3 ■ 








-4.7 


-4.4 


-4.0 


-2.4 


-4.6 


-3.7 


-2.6 


10 3 • # 2 








1.4 


1.2 


1.2 


1.0 


2.3 


-0.2 


0.8 




768.3 


0.5 


100 


811 


830 


831 




802 


1260 




M A 


1260 


30 


300 


1331 


1376 


1609 




1610 


2010 




fv 


0.20 


(*) 


0.02 


0.18 


0.17 


0.17 




0.18 


0.15 




9v 


0.090 


(*) 


0.009 


0.081 


0.079 


0.079 




0.080 


0.076 




fA 


0.097 


0.022(*) 


0.022 


0.083 


0.080 


0.068 




0.072 


0.084 




M s 


983.3 


2.6 


200 


617 


620 


709 


989 


657 


643 


760 


Cm 








20 


18 


20 


24 


25 


16 


6 


Cd 


34 


(*) 


10 


21 


21 


18 


23 


19 


26 


27 


X 








0.052 


0.063 


0.057 


0.089 


0.035 


0.1 


0.2 


9a 








0.61 


0.62 


0.62 


1.0 


0.66 


0.79 


1.0 


Mq 








265 


263 


246 


199 


204 


262 


282 


9 








0.0 


0.0 


0.25 


0.58 


0.5 


0.0 


0.0 



Table 1: Experimental values and predictions of the ENJL model for the various low-energy 
parameters discussed in the text. All dimensionful quantities are in MeV. The difference 
between the predictions is explained in the text. The numerical error in |13| for H2 has 
been corrected. All masses are determined from the low-energy expansion, not the pole 
position of the 2-point functions. 
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9 Two-point functions 



This section is a discussion of the results in Refs. |53|, [54]] about two-point functions. These 



two-point functions were studied before in but there they were discussed as quark form 
factors. What is new here is that the explicit dependence on the regularization scheme has 
been put into two arbitrary functions, namely, Uy 1 + Ily and LT^ (see this section below 
for definitions). This also shows that these results are valid in a class of models where the 
one-loop result can be expanded in a heat-kernel expansion using the same basic quantities 
E and as used here. This includes the ENJL model with low-energy gluons described 
by background expectation values. 

The two-point functions are of course important quantities and have played historically 
an important role in understanding the high-energy behaviour of the strong interaction ||100| , 



1 6|1 . In addition some of the consequences for the mesonic sector were also valid in the 
low-energy expansion of the ENJL model as discussed in the previous section. Here we 
would like to study the two-point functions directly in the ENJL model to all orders in the 
current quark masses and momenta. This method was developed for the chiral limit case 
53] and then extended to include nonzero quark masses in |54|. The discussion here 



m 



follows the latter reference closely. 

9.1 Definition of the two-point functions 

We shall discuss two-point functions of the vector, axial-vector, scalar and pseudoscalar 
quark currents with the following definitions, 

V*{x) = UxhMx), ( 209 ) 

JTfa) = qi(x) Wqj (x) , (210) 

S ij (x) = -qi(x)qj(x), (211) 

PV{x) = qi(x)i^ qj (x), . (212) 

The indices i,j are flavour indices and run over u, d, s. The two-point functions themselves 
are defined as 

KMm = ^/d 4 xe^<0|T(^(a;)K fe/ (0))|0>, (213) 

KM^ki = z/d 4 a;e^<0|T(4(a;)^(0))|0>, (214) 

n£(?)tf« = i J d 4 xe"^ < 0\T (v» (x)S M (0)) |0 > , (215) 

K^ki = * J dWf* < 0\T (^(x)P fc '(0)) |0 > , (216) 

Tl s (q) ljH = i J d A xe^ x < 0\T (S i \x)S H (0)) |0 > , (217) 

U p {q) ijH = ijd 4 xe iq - x < 0\T (P ij {x)P H (0)) |0 > . (218) 
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In the leading order in the number of colours these are all proportional to Sij k i = SaSjk) with 
Su the Kronecker delta. Using Lorentz-invariance these functions can then be expressed as 
follows 



KMm = {(^-?V) n v(^ 2 )y + ^ n y(^ 2 k}^> ( 219 ) 

KMvti = {(q»qu-q 2 g»u)n%\Q% j + q„q u nf(Q 2 ) ij }5 ijkl , (220) 

K%q)i 3 ki = q^f{Q 2 )ij8 ijkU (221) 

nj(g) yw = iq^{Q%8 ijkU (222) 

U s {q) ijkl = ILs(Q%6 ijkl , (223) 

U p (q) m = U P (Q%5 tjkl . (224) 



Here Q 2 = —q 2 . We shall discuss the Weinberg Sum Rules and numerical results for the 
two-point functions only in the Euclidean domain, i.e. Q 2 positive. Using Bose sym- 
metry on the definitions of the two-point functions it follows that Uy\Q 2 )ij, ILy (Q 2 )ij, 
l4 0) (Q%, U^(Q 2 ) l3 , Ii s {Q 2 ) l3 and U M (Q 2 ) ij are all symmetric in the flavour indices i 
and j. The remaining ones need the Ward-identities to prove their flavour structure. From 
the identities in the appendix [A] it follows that Hg(Q 2 )ij is also symmetric in while 
^¥(Q 2 )ij i s ant i- symmetric. 

9.2 Lowest order results in Chiral Perturbation Theory 

From Chiral Perturbation Theory to order p 4 in the expansion we obtain the following low 
energy results for the two-point functions. The orders mentioned behind are the orders in 



Chiral Perturbation Theory that are neglected. 

n$\Q% = -4(2ffi + Lio) + 0(p 6 ) , (225) 

4 0) (Q%- = 0(p e ), (226) 

n?(Q%' = - 4(21^ - L 10 ) + 0(p e ) , (227) 

Kf(Q% = 2/ 3 (- f l^-i,) +0 (p 6 ), (228) 

Tlf(Q% = 0(p 6 ), (229) 
IB f 2 

n2W% = ^T^ + (^ 6 )' (230) 

UsiQ 2 )^ = 8B 2 (2L 8 + H 2 ) + 0(p 6 ), (231) 

Kp(Q 2 h = 2 2°!% + 8B o(~ 2L s + H 2 ) + 0{p*) . (232) 



With rriij the mass of the lightest pseudoscalar meson with flavour structure ij. These are 
obtained in the leading 1/N C approximation so loop-effects are not needed. Notice that 
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these expressions are valid to chiral order p A . From a term of the form ^{D^xD^x^} there 
are contributions of order (m, — rrij) 2 /Q 2 to the vector two-point function Uy\Q 2 )ij and 



of order (m^ — rrij) to the mixed scalar vector function Ug^Q 2 ^ 

The functions n^, Hp and lip get their leading behaviour from the pseudoscalar 
Goldstone pole. In addition and contain a kinematical pole at Q 2 = 0. The 
residue of the physical pole is proportional to the decay constant for the relevant 
meson, (for the ud ones, f u d ~ ~ 92.5 MeV). In xPT, the constant B is related to 
the vacuum expectation value in the chiral limit. In the large N c limit and away from the 
chiral limit there are corrections due to the terms proportional to combination of 0(p A ) 



couplings 2L 8 + H 2 |p9 



< 0| : : |0 >\y= uAs = -f*B (l + 0(p A )) . (233) 

The vacuum expectation value here, < 0| : : |0 >, is the one used in %PT in the chiral 
limit and fo is the pseudoscalar meson decay constant in the chiral limit. The constants 
L 8 , L w , Hi and H 2 are coupling constants of the 0(p A ) effective chiral Lagrangian in the 
notation of Gasser and Leutwyler [ffSP , section [|. The constants L 8 and L 10 are known 
from the comparison between %PT and low energy hadron phenomenology. At the scale 
of the p meson mass they are L 8 = (0.9 ± 0.3) x 10 -3 and L 10 = (—5.5 ± 0.7) x 10~ 3 . The 
high energy constants Hi and Hi correspond to couplings which involve external source 
fields only and therefore can only be extracted from experiment given a prescription. 



9.3 The method and Ward identities 

The method used here is identical to the one used in The full two-point functions are 
the sum of diagrams like those in figure []a. The one-loop two-point functions are those 
obtained by the graph in figure |]b. Using a recursion formula that relates the n-loop graph 
to a product of the one-loop and the (n-l)-loop graph and the relevant combination of 
kinematic factors and Gy and G$ the whole class of graphs can be easily summed. Some 
care must be taken in the case where different two-point functions can mix so a matrix 
inversion is necessary (see ref. |53| ). 

The two-point functions defined above satisfy the following Ward identities. (We sup- 
press the argument Q 2 for brevity.) 

K-m,)n^., (234) 
(mi - rrij) U Sij + faqi) - (qflj) , (235) 
K + m^n^., (236) 
(mi + rrij) Tl Pij + (q^) + (q^j) . (237) 

These are derived in the appendix 0. From these the flavour symmetry of the mixed 
two-point functions can be derived from the vector ones. 



-W J7 = 
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(a) (b) 



Figure 6: The graphs contributing to the two point-functions in the large N c limit, a) The 
class of all strings of constituent quark loops. The four-fermion vertices are either £njl 
or jC^jl in eqs. (|14|) and ([15]). The crosses at both ends are the insertion of the external 
sources, b) The one-loop case. 



The one-loop expressions, which we shall denote by II and use further the same con- 
ventions as given for the full ones above are given in appendix [FJ. They satisfy the same 
identities but with the current quark masses m; replaced by the constituent ones, Mj. 
In addition to these, there are two more relations that follow in general if the one-loop 
part can be described by a heat-kernel expansion in terms of the quantities E and of 
appendix |A|. These identities are (with the flavour subscript ij and argument suppressed) 

n«+4 0) = n« + ni 0) , (238) 

U S + Q 2 U^ = Tip + Q 2 T[f. (239) 



Let us show in the most simple case|66| how there are general consequences of this 



approach. We will derive here the relation between the scalar mass and the constituent 
quark mass in the chiral limit. The set of diagrams that contributes is drawn in Fig. |6]a. 
The series can be rewritten as a geometric series and can be easily summed in terms of the 
one- loop 2-point function II5. The full result for the scalar-scalar two-point function (we 



only treat the case with equal masses here, see [54] for the general case) is: 



n 5 = : Us = ■ (240) 
1 - g s IL s 

The resummation has generated a pole that corresponds to a scalar particle. Can we say 
more already at this level? 

We can in fact. The Ward identities for the one loop functions become: 

n s = u P - q mf, (241) 

Tt _ <? 4 ^(0) (qq) (9d? , 

Up ~ mf A ~ M~ Q • (242) 
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( |241|) is a consequence of using the heat kernel for the one-loop functions and (|242 ) is a 
direct consequence of the symmetry. Using these two relations we can rewrite 



2tt(°) 



The first two terms vanish due to the gap equation so this two-point function has a pole 
at twice the constituent mass. Let us now derive the general cases. 

9.4 The transverse vector sector 

We use here the abbreviations (|1^). The full resummed transverse vector two-point function 



is then 

4\ = -< . (244) 
1 + g v n Vij 

This can be simply written in a form resembling the one in the complete VMD limit with 
couplings fs, fv an d My depending on Q 2 and flavour and defined by 

rti^tn^ - 2 fs(Q 2 )ij | 2 fv(Q 2 )ij M v(Q 2 )ij ( 0A r\ 

Ey(gk " ~^^ + M 2 {Q 2 ) l3 + Q 2 ' (245) 

lfl(Q 2 h = ~ Q2U ^ 2 } j , (246) 



9vK\Q 2 : 



'j 



2fKM(Q% = ^J.,^ , (247) 
2f v {Q% 3 = T^ +1) {Q 2 ) l3 . (248) 

Where we have used the fact that (see appendices [A] and |B|) 11^ 1 = Ily + li^y has no 

pole at Q 2 = 0. There is a correction here (in Uy ) due to the mixing with the scalar 
sector, which is allowed by the presence of explicit breaking of the vector symmetry (see 
the scalar mixed sector subsection |9.7[). For the diagonal case, this is defined as m ; = m,- 



or Mi = Mj, Uy vanishes and the formulas above simplify very much. 

The pole mass of the vector corresponds to the pole in this two point function or to 
the solution of Re (Q 2 + M v (Q 2 )ij) = 0. Alternatively, one can define the VMD values for 
the vector parameters (fy and My) as the best parameters of a linear fit of the inverse of 
Tly (Q 2 )ij — 2fs(Q 2 )ij/Q 2 . These definitions have the advantage that they are also valid for 
the Euclidean region (Q 2 > 0) where the vector cannot decay into two constituent quarks. 



See sections on numerical applications |9]9] and Vector- Meson-Dominance |10.3| for further 
comments. 
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9.5 The transverse axial-vector sector 



The transverse axial-vector two-point function derivation is also identical to the one in ref. 

= (249) 

A ii -> . rf(i) ' 1 yj 

1 + gv^-Aij 

Using the identity (J238|) it can be seen that this has a pole at Q 2 =0 because 11^ has it. 



As can be seen from the explicit expression and is proved in general in appendix 0, the 
the pole at Q 2 = in a simple fashion. 



combination Ily' ) + ILy is regular at Q 2 going to zero. This again allows us to separate 



u^ QX = « + , (250) 



Q 2 M 2 A (Qi) tJ + Q 2 

fijiQ 2 ) = 9A(Q%f*(Q 2 ), (251) 
2jJ(Q 2 ) = -Q^^k , (252) 

^(g 2 )ii) _1 = l-^n^^, (253) 
iV A 2 

2/ 2 (Q 2 ),,M 2 (Q%. = — -|_2L^ A (Q 2 ) i . , (254) 

2/l(Q 2 )i,- = g 2 A {Q\^v +1 \Q 2 h ■ (255) 

There is a correction here (in 11^) due to the mixing with the pseudo-scalar sector due to 
the presence of both spontaneous and explicit breaking of the axial- vector symmetry (see 
the pseudo-scalar mixed sector subsection). For further discussion of these expressions and 
the ones in the previous section we refer to the subsection pT8| on Weinberg Sum Rules. 



9.6 The pseudo-scalar mixed sector 



The results obtained in |53|, [54] are, with the summed functions given in terms of the func- 
tion A P (Q 2 ) and the one loop two-point functions (with flavour subscripts ij suppressed), 



n? } (Q 2 ) 

n P (Q 2 ) 
Ap(g 2 ) 



A P (Q 2 

_J_pM 

A P (Q 2 ) 
1 



A P (Q 2 



1 - g s Tlp(Q 2 ))Tif(Q 2 ) + gs(Up(Q 2 )f 

np(g 2 ), 

1 - gyllf (Q 2 ))TT P (Q 2 ) +g v (Jlp(Q 2 )f 



l-g v uf(Q 2 ))(l-g s U P (Q 



gsgv 



(256) 
(257) 
(258) 
(259) 



Using the identities for the one-loop case it can be shown that the resummed ones satisfy 
the Ward identities of appendix |A| with the current quark masses. To show this it is also 
necessary to use the Schwinger- Dyson equation for the constituent quark masses in eq. ( 
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In order to rewrite this in terms of a nicer notation we first express Ap(Q 2 )ij in a 
different form using the identities for the one-loop two-point functions. 



A P (Q% 
with m 2 AQ 2 ) 



Mi + Mj 



{ml(Q 2 ) + Q' 



(mi + mj) 



■M , 



g s g A (Q 2 )Tl P \Q 2 



(260) 
(261) 



Inserting the definition of fjj{Q 2 ) and 1/gs = 
Mann-Oakes-Renner (GMOR) relation for the pion mass 



(liQi) I (Mi — rrii) we recover the Gell- 
when eq. (|261|) is expanded 



in powers of mj. For further discussion on corrections to the GMOR relation in this model 
we refer to the section on numerical applications |9.9| . Formula Q261| ) gives the expression for 
the pole due to the lightest pseudoscalar mesons in the presence of explicit chiral symmetry 
breaking. 

This then allows us to rewrite the full two-point functions in a very simple fashion: 



np*(Q 2 



Mi + Mj 



™UQ 2 ) 



Q 2 Q 2 



9s m 2 (Q 2 )+Q 2 ' 
1_ (Mj + Mj) 2 1 1 
9s 2MQ 2 ) g 2 s m 2 3 (Q 2 )+Q 2 ' 



(262) 
(263) 
(264) 



Here we want to point out that the two-point functions Up and lip suffer from the 
same ambiguity (via its dependence on g$) as the quark-antiquark one point-function (see 
discussion at the end of section ^) when compared with the %PT results. 



9.7 The scalar mixed sector 

This can be done in the same way as in the previous subsection with the result (with 
flavour subscripts ij suppressed)! 



4 0) (Q 2 ) 
nf(Q 2 ) 

n 5 (Q 2 ) 

A S (Q 2 ) 



1 



As(Q 2 ) 

1 nf(Q 2 



1 - g s U s (Q 2 ))Il^\Q 2 ) + g s {Tlf{Q 2 )) 2 



A 5 (Q 2 
1 



(1 -g v U^(Q 2 ))U s (Q 2 ) +g v (uf (Q 2 )) 



1 - g v nf(Q 2 )) (l - g s Tl s (Q 2 )) - 9s9v(^(Q 



(265) 
(266) 
(267) 
(268) 



To rewrite this in a simple fashion we would again like to expand As in a simple pole like 
fashion. Using the identities for the one-loop two-point functions this can almost be done, 
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we obtain 



A S (Q 2 )« = Sj 0^f {(Mi + Mjf + ^(Q 2 )«mJ(Q 2 ) + Q 2 ) 

It can be seen that in the diagonal case a simple expression for the scalar meson pole can 
be found, 

(Mi + Mj) 2 + g A (-M 2 s ) u ml(-M 2 s ) . (270) 



M 2 S (-M 2 ) 



mi=m,j 

2\ 



The expression for the scalar two-point function Tls(Q ) is in this case 



So in the diagonal case a simple relation between the scalar mass, the constituent masses 

r(0) 



and the pseudoscalar mass remains valid to all orders in the masses. In this case 



Uf = 0. 



For the off-diagonal case, i.e. rrii ^ mj, the corresponding expressions for nj?, 
and Us can be obtained from eqs. ( |265| )-( [2"68| ) and the explicit II functions in appendix 
[FJ. There is a small shift in the pole compared to eq. ( |270|) for the case mj ^ rrtj. From 



appendix [B|, in eq. ( |B.4j) , it can be seen that li^y itself has a zero close to a value of 

Q 2 = M 2 S of eq. (E7Cp - In addition TT^ is suppressed by (Mj - Mjf /Q 2 . Therefore the 
value of the pole in the off-diagonal case is not too far from that in eq. ( |270| ). 

Here we want to point out that (as in the mixed pseudoscalar sector) the two-point 
functions Tly\ Iljf and lis suffer from the same ambiguity (via its dependence on g$) 
as the quark-antiquark one point-function (see discussion at the end of section |2|) when 
compared with the xPT results. 



9.8 Weinberg Sum Rules 

The Weinberg Sum Rules are general restrictions on the short-distance behaviour of various 
two-point functions [|lQQf| . They were first discussed within Q CD in ref. |[L06|1 . A low-energy 
model of QCD should have a behaviour at intermediate energies that matches on reasonably 
well with the QCD behaviour. The general behaviour should be (II^r = ITy — IT4.) 

FirstWSR , (272) 
SecondWSR , (273) 
ThirdWSR . (274) 



hm Q 2 niTW 

Q — >oo N 

hm (g 4 n«(Q 2 
hm (g 4 iffi(Q 2 
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Let us review first the QCD behaviour of these Sum Rules. In the large N c limit the three 
WSRs are theorems of QCD in the chiral limit (i.e., M. — > 0). The first WSR is still fulfilled 
in the large N c limit with non-vanishing current quark masses. However the second and 
the third ones are violated as follows 
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lim (Q 4 IlS(Q 2 



-JmjQ 4 nS(Q 2 ) 
2 (rriifaqj) + rrij(qiqi)) 



(275) 



As shown in [53] the class of ENJL-like models does satisfy the three WSRs in the 
chiral limit. We shall now check how well this does in the case of explicit breaking of chiral 
symmetry. 

The high-energy behaviour of the two-point functions ITy^ needed for the three WSRs 



can be easily obtained from the expressions in sections 9.4, 9.5, |9.6| and 9.7. The first 



and second WSRs are satisfied in these ENJL-like models even with non-vanishing and all 
different current quark-masses. The high energy behaviour (Q A ) of these models is thus 
too strongly suppressed for n^j^Q 2 ) to reproduce the QCD behaviour in the second WSR. 
The third one is violated as in QCD and one has 



lim (Q 4 nS(Q 2 



9s 



{rriiMj + rrijMi) 



(276) 



Let us now see what relations between low-energy hadronic couplings do these Sum 
Rules imply for this ENJL cut-off model. In the equal mass sector, rrii = rrij ^ 0, one has 



f 2 M 2 = f\M\ + fl 

fyMy = f\M\ 



Remember that in QCD one has in this case 



fyMy 



f A M\ + fl 



f v M v = f A M\ + miti . 



(277) 
(278) 



(279) 
(280) 



In the off-diagonal case, mj ^ rrij, the situation becomes a lot more complicated. However, 
since the off-diagonal part is suppressed by (Mi — Mj) 2 /Q 2 one does not expect qualitatively 
different results. 



9.9 Some numerical results 



As can be seen from the explicit formulas the change with respect to ref. [|53[] is in most 
cases a (small) shift in the two-point function mass pole positions. Therefore we do not 
plot too many of the two-point functions. As numerical input we use for G$, Gy and A x 

ITSfl. These are A v = 1.160 GeV and Ga = 1.216. The value 



the values from fit 1 in ref, 
°f 9a(Q 2 = 0) there was 0.61. This is Gy 
the value of the quark mass for m = m u - 



- 1.263. For the current quark masses we use 
irid that reproduces the physical neutral pion 
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m1=m2=0 
m1=m2=3.2 
m1=m2=83 




-0.5 



Figure 7: The inverse of the transverse vector two-point function for equal quark masses 
in the chiral limit, i.e. M. — > 0; for the p meson, i.e. m\ = vri2 = 3.2 MeV and for the <f> 
meson, i.e. mi = rri2 = 83 MeV. The units of q 2 are GeV 2 



and kaon masses. With the other parameters as fixed above this is m = 3.2 MeV and 
m s /m = 26. 

As an example we have plotted the inverse of the transverse vector two-point function 
in eq. ( p45| ) in figure |7] for the values of Gs and A x corresponding above mentioned. The 
full curve is the result in the chiral limit (M. — > 0) and the dashed is the result with 
rrii = rrij = m the value above. The reason we have plotted the inverse will become clear 



in section |10.3| . We also show the inverse for rrii = rrij = m s the value above in the short- 



dashed curve. To show the result for unequal quark masses we have plotted in figure || 
the transverse vector two-point function itself for the chiral limit case and for the us case 
with m s and m above. Notice that the two-point function now has a kinematical pole at 
q 2 = 0. 

We have also plotted in figure |^ for the parameters quoted above the dependence of the 
pion mass on Q 2 . Since /^m?- is a constant, see eq. ( |261| ) this is also the Q 2 dependence 
of the inverse of the decay constant squared. 

Let us make some comments on the corrections we find to the GMOR relation ( [26 1|) 



in this model. The corrections to the GMOR relation ||108|| can be calculated here in an 
analogous expansion to the one in x?T. Then the GMOR relation can be written as follows 
(for the diagonal flavour case, i.e. i = j) 



(281) 
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-0.15 - | 

-0.2 I ' 1 1 1 h 1 1 

-0.5 -0.4 -0.3 -0.2 -0.1 0.1 0.2 

q A 2 

Figure 8: The transverse vector-two-point function for the chiral limit and for unequal 
quark masses, m 1 = m and m 2 = m s . Note the kinematical pole at q 2 = 0. The units of 
q 2 are GeV 2 . 




0.015 



q A 2 (GeV A 2) 



0.15 



Figure 9: The running pseudoscalar mass squared, m 2 A—q 2 ), as a function of q 2 for = 



rrij = 3.2 MeV. 
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Here we have included all the chiral corrections to the quark condensate, to the pion 
mass and to the pion decay constant in their respective values. Then the remaining is a 
correction to the GMOR relation. We have also calculated this correction in this model 
and it turns out to be 

1 - 4^ (2L 8 - H 2 ) + 0{f) = A (l - ^) . (282) 

Notice that the r.h.s. contains all the orders in the ;^PT expansion in the large N c limit. A 
is the reducing factor discussed in eq. ( |312| ). Numerically, this correction is around 0.2 % 
for pions and 3 % for kaons and disagrees with the one found in QCD Sum Rules, see the 
book by Narison in and [FT . The expression in (|282|) is the one consistent with the use 



of Ward identities to sum the infinite string of constituent quark bubbles. Numerically we 
get 2L 8 — if 2 — 0.2 ■ 10~ 3 for the input parameters above. This agrees with the one found 
at the one- loop level, Sect. 0. 

9.10 Inclusion of gluonic effects 



The inclusion of extra gluonic effects like in || [L3|, |53| can be done simply by replacing 
the one-loop functions of appendix IB] by the ones including the effects of (asG 2 ). the 



expressions for the chiral case can be found in ||53|| , these can also be used for the case of 
equal but nonzero current quark masses. In the general case the expressions needed can 
be found in the QCD sum rules literature, see e.g. ||]. 

The effects are in general rather small because the main behaviour is produced by the 
resummation and not so much by the slow variation with Q 2 of the one-loop functions. An 
example was shown in Fig. 4 of fl5"3f . 



10 Some results on three-point functions and Meson 
Dominance 

In general the same procedure as used in the previous section can be extended to three-point 
functions. In fact most of them were used in the calculation of the Bk factor [^, 60| . Here 



we will discuss one example extensively and a second one in a short form. Some comments 
about meson dominance in the three-point functions are also given. The discussion closely 
follows Ref. 



10.1 VPP with the use of the Ward identities 

In this subsection we calculate the Vector Pseudoscalar Pseudoscalar (VPP) three-point 
function to all orders in %PT using the same type of methods as those used for the two-point 
functions. The three-point function we calculate is the following 

n^ PP (pi,p 2 ) = i 2 J d 4 x J dV (pi ' x+P2 ' y) (0|T (v; j (0)P kl (x)P mn (y)) |0) . (283) 
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Where i,j, k, I, m and n are flavour indices. In the limit of large N c the flavour structure is 
limited because of Zweig's rule (this flavour structure is general for any three-point function 
of three quark currents), 



n ^ PP (Pl ; P2) = ^(jPl,Pz)ikm3ilSknS m j + ^ (pi, P2)ikm$inhj$ml ■ (284) 

Bose symmetry requires that 

n^(pi,£>2)ifcm = n.~(p 2 ,Pl)imk • (285) 

The three-point function il^ pp (pi,p2) can then be simply calculated by only taking one 
particular flavour combination. Finally we can use Lorentz-invariance to rewrite 

nj(pi,p 2 )ifem = Pi,tllt km (pl,pt,q 2 ) +p 2 „Ilf km (pl,pt,q 2 ) , (286) 

where we have defined q = pi + p 2 . 

We shall limit ourselves to the vector diagonal case, i.e. rrii = rrij. In the vec- 
tor off-diagonal case there will also be non-trivial mixings with the scalar-pseudoscalar- 
pseudoscalar three-point function. Here a relatively simple Ward identity for this three- 
point function can be derived from d^V^ = and the equal-time commutation relations. 
It is 

q ,1 n^( Pl , P2 ) lkm = -Upi-pDki + U P {-p 2 ) mk . (287) 

So the Ward identity relates the three-point function to a combination of two-point func- 
tions. This determines one of the two functions IT 4 , U B in terms of the other. The Ward 
identity gives, for instance, the following constraint (for p\ = p\ and i = m) 

nf fc ,(p 2 ,P 2 ,g 2 ) = -ni(p 2 ,p 2 ,g 2 ) . (288) 



The type of graphs that need to be summed are depicted in figure [L0| Each of the 
three tails here is the diagram in figure §a with the same explanation as there. We have 
there depicted one particular flavour combination. This is the one that corresponds to the 
function 11+ given above. The i, k, m written above the lines are the flavours of each line. 

All graphs are formed by having the tails summed over 0, 1, 2, • • •, oo loops connected 
by four-fermion couplings. These then couple to the one-loop three-point function (or 
vertex) II , with various possibilities for the insertion in the three-point vertex. These 



possibilities for the 7-matrices are written in figure |10| inside the main loop. 

In this figure the left-hand side depicts the insertion of the current V^(0) and Tail I is 
the connection to this current. On the end connecting to the one-loop three-point function 
it is only nonzero for another vector insertion since in the diagonal case we consider, the 
mixed vector-scalar two-point function vanishes. Its expression is given by 



!!,»■ + ^ n£(-g 2 ) mi * ( 289 ) 
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Tail II 



Tail I 



7/, 




m 




Pi 



m 



m 

Tail III 



275 

P2 



Figure 10: The graphs that need to be summed in the large N c limit for the Vector- 
Pseudoscalar-Pseudoscalar three-point function. See text for explanation. 



Here the first term comes from where the external current directly connects to the one-loop 
three-point function and the second term is with the two-point function in between. The 
sum of both is 

A similar discussion can be done for Tail II and Tail III. First we have the insertion 
of the current P kl (x) at the external end. On the end connecting to the one-loop three- 
point function we can have 275 or an axial-vector insertion since the mixed axial-vector- 
pseudoscalar two-point function is nonzero. The 275 insertion tail is : 



1 + 



np(-Pi)fci 



(M fc + Mi 



29sf 2 ki (-pi)(™U-pi)-pl) ' 
For the connection with the axial-vector insertion it is instead 



(291) 



N c Kl 



ki 



(M k + Mi 



2tiM>g s (mU-pl)-pj) 



(292) 
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The combination f v M v is here flavour and p\ independent, it is the combination in eq. 
( j247|) with Hy\Q 2 )ij = since we are in the diagonal flavour case. The way both these 
types of insertions can appear due to the tail are how within this formulation the mixing of 
pseudoscalar and axial- vector degrees comes about. These will be described by factors of g\ 
(see below). Tail III is identical to Tail II with the substitutions p\ — > P2 and i, k — > k, m. 
The full expression for 11+ is 

n + "( Pl , P2 ) = {<r + z ^ Kn " 7 "'(-9 2 )™} 

+n™ (Pl , K) (i + ^n P (-rf) t> ) (^rfn?(-riu) (293) 

YPA VAP VAA 

Where the one-loop three-point functions II v , II and II va are the one fermion-loop 
result for 



Ul u PA ( Pl ,p 2 ) = i 2 j d 4 x J dV (pi " +P2 ^(0|T(v; m (0)P fci (x)^ fc (y)) |0>, 

(294) 

nl„ AP ( Pl ,p 2 ) = i 2 j d 4 x J d 4 ye^ + ^(0\T(v^(0)A k ;(x)P mk (y)) |0> , 

(295) 

n^(pi,p 2 ) = i 2 J d 4 x y dV^+^OlT (^(0)^(x)A- fc (2/)) |o) . 

(296) 



To obtain the full expression in eq. ( |293| ) it now remains to calculate these VPP, VAP, 



VPA and VAA one-loop three-point functions (or vertices). The axial- vector ones always 
come multiplied with the relevant momentum. So we always have the scalar products 
Pi ■ A ht (x) and p 2 ■ A mk (y). That means that using the Ward identities we can relate the 
VAA, VAP, VPA to the VPP one plus possibly two-point function terms resulting from 
equal time commutators. 

These Ward identities are (remember we assume Mj = Mj here). 

irfn%?(pi,P2) = -(M k + M t )Tl V ^ A (p 1 ,p 2 ) 

+*njl(-g) mi - iU A a (-p 2 ) mk ; (297) 

i P »Tilf(pi,p 2 ) = -(M fe + M,)n+( Pl , P2 ) 

+iTlp^-p 2 ) mk . (298) 
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The other needed ones can be derived from this using Bose-symmetry. Notice that there 
is no contribution here from the flavour chiral anomaly. 

We can now use these identities to obtain the final result for the three-point function 
we want. The terms which after the use of the one- loop identities above are proportional 
to VPP can be combined into a simple form using <7aG° 2 )- The result is (we have Mj = Mj 
and j = m in this flavour configuration) . 



n + ^(pi,p 2 ) 



(Mj + M k f \ f g ^M 2 (-q 2 ) mi -q^ 
MfU-Pi)fL(-Pl)) V M 2 (-q 2 ) mi -q 2 



Jl-9A(-pl)»Hl-9A{-tiU) {{p2 . ^ _ {pi . q)p2u}T[ U { _ q z U (299) 



(Mi + Mj 



k) 

gA{-p\)ki (1 - 9A(-Pl)mk) jtM, 2 \ 

M . + Mk P^P \ Pi)* 

g A {-pl)mk {l- 9A(-Pi)ki) jjM, 2 . 



Mi + M k 

This result satisfies the Ward identity fl287p if the one-loop function one satisfies the 
same one with the one-loop functions. This provides a rather non-trivial check on the 
result 



It now only remains to calculate the one-loop form factor U (pi,p 2 ). We give its 
expression in appendix |C|. At this point we can see in eq. ( |299| ) how far regularization 
ambiguities affect the result. We first have to define the two-point functions. Here all 
ambiguities are restricted to two bare functions (see section [5] for details). This three-point 
function adds one more in general, the three-propagator function 7 3 (pf,p|, q 2 ) (see explicit 
expression in appendix 0). Of course, this one-loop form factor EL (pi,p2), satisfies all the 
identities eqs. ( |283| ) to Q288Q as well. We refer to section [L0.3.2| for the definition of the 



physical vector form factor after reducing this VPP three-point function. We shall also 
discuss there the VMD limit in this form factor and give some numerics. 

The same three-point function can be calculated in Chiral Perturbation Theory. The 
result is 

- «i - - Pi) fa - ^ r i 1 + + om ) ' (300) 



Pulling out the pion poles (see section |10.3.2| for technical details) and taking the low- 



energy limit and the value of L 9 in this class of models our full result in eq. (|299|) reduces 
to this, providing one more non-trivial check. 
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10.2 PVV with a discussion about its Ward identity 

In this subsection we calculate the Pseudoscalar Vector Vector (PVV) three-point function 
to all orders in ;\PT with the same method as the one used before. 

n^ y (pi,p 2 ) = i 2 J d 4 x J dV^+^OlT (P^(0)V* l (x)Vr (y)) |0> . (301) 

The pseudoscalar can couple with at the one-loop end, both an axial-vector and pseu- 
doscalar two-point function. These have the same form as equations ( |291| ) and (292) in 



the previous section with p\ —> q. Summing up the three tails the total result is then [54 



n+^( P i,p 2 ) = 



X 



J^AT" ■ ™>)\" ^^Af 11 ™; (302) 



x 




with the one-loop result for the following three-point function 



nfvv( pi , p2 ) = i*J d 4 x J d A ye^ + ^y\0\T(A l p m (0)V^\x)V u mk (y)) |0) . 

(303) 

The main new part here is that at the one-loop level we now have to include the 
anomalous part of the Ward identities. There has been in fact quite some confusion whether 
this can be done consistently, see Sect. [7| Here we want to apply that method to the PVV 
three-point function to all orders in external momenta and quark masses. The prescription 
is essentially to use the anomalous QCD Ward identities for the axial current consistently. 
We shall use the scheme where vector currents are conserved ||93|| . When we use the one- 
loop anomalous Ward identity to reduce the right-hand side of the pseudoscalar to a part 
with only pseudoscalar couplings to the one-loop vertex, we obtain a local chiral invariant 
result plus an extra part where the tail couples directly to the external vector sources 
v kl (x)v™ n (y) . This extra part is of order p 4 and is the subtraction the anomalous Ward 
identity imposes to obtain the correct QCD flavour anomaly. 

The full result in terms of the one-loop n^, three-point function is given by 



n+(pi,p 2 ) = n^(pi,p 2 ) 



(mU-pI)u-p?) (mU-pD* 



'l + ^n P (-g)-^2MT^ 

AT A 2 r\ 1) AT A 2 1 f 



N C A Z X y 11 N C A X 



+ n (pi,p 2 ) 



8tt 2 G 

Pi=p 2 2 =ql=0 NcA X 



v 2MiUf \-q) . (304) 
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Where the one-constituent quark loop function U is given by 



with F(pl,plq 2 ) = l + I 3 (pi,plq 2 )-I 3 (0,0,0) (305) 

where the form factor l3{p 2 ,P2,q 2 ) is the one given in appendix y and which appeared 
before in the study of the VPP three-point function in section ( |10.1|) . This form factor 
coincides with the one found in the context of constituent quark-models (see for instance 
109|| ) when the cut-off A x is sent to oo. Here, this is a physical scale of the order of 



the spontaneous symmetry breaking scale and therefore we have to keep it finite. The 
anomalous Ward identities are telling us that terms which are of chiral counting different 
to 0(p 4 ) have to be local chiral invariant but they do not fix the regularization for those 
terms. We therefore use here consistently the same regularization for them as in the non- 
anomalous sector. At 0(p 4 ) the chiral anomaly also uniquely fixes the one-loop constituent 
chiral quark anomalous form factor to be the one in eq. ( |305|) when p\ = p\ = q 2 = ||110 1 . 



Here we have used the anomalous Ward identity in eq. ( [30 2| ). A naive use of the two- 
point functions and Ward identities would have led only to the first term in the sum in eq. 
( |304|) . The second term is the result of enforcing the validity of the QCD flavour anomaly. 
Substituting the results on the two-point functions in section |^ we can write down the 
following explicit expression 



167r \9sfu(-q 2 ) {mli-q 2 " 



9A(-q 2 ) 



i F( J „2 2, M 2 (~plhM 2 (-plh 
1 - r (Pl)P2) Q )- 



(M 2 {-p\) u -pl) (m 2 {-pI) u -pI 



(306) 



We refer to section |10.3.3| for the definition of the physical anomalous 7r°7*7* form factor 
after reducing this PVV three-point function. We shall also discuss there on the VMD 
limit in this process and give some numerics. 

10.3 Meson-Dominance 

We already saw that in the low energy limit we had a lot of relations that were equivalent 
to various meson dominance relations. Here we discuss the extension of those relations to 
the all-order case. 

10.3.1 Two-point functions 

Here we shall discuss the vector case, the axial- vector case is similar. The transverse vector 
two-point function in eq. ( |245|) reduces in the diagonal case, rrii = rrij (the off-diagonal 
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case can be done analogously) to the following simple expression 



nPhrt = m(-, 2 ) M f { l[ 2 fl q2 (307) 

N A 2 

with 2/ 2 (-g 2 )M 2 (-g 2 ) = _£-2L (308) 

o7T (jy 

and 2/ 2 (-g 2 ) = ^\-q 2 ). (309) 

In the complete VMD limit this two-point function has the same form but with fy and My 
constants. Let us see how complete VMD works in this model. For that, we shall study 
the inverse of Uy\—q 2 ), which in the complete VMD limit is a straight line. This function 
was plotted in section |9]9| in figure [7]. There we can see that U ( y\ — q 2 ) in this model is very 
near to reproducing the complete VMD linear form. Moreover, we can perform a linear fit 
to the inverse of Tly\—q 2 ) to obtain the best VMD values for the fy and My parameters. 
These parameters are in this way meaningfully defined in the Euclidean region — q 2 > 
where the model is far from the two constituent quark threshold. Doing this type of fit for 
the values of the input parameters A x , Gy, Gs discussed in section |9]9] leads to My ~ 0.644 
GeV for the vector mass in the chiral limit (remember that we are always in the large N c 
limit) and fy ~ 0.17 for the decay constant. For current quark masses values discussed 



also in section |9.9| , we obtain for the p meson flavour configuration M p ~ 0.655 GeV and 
f p ~ 0.17 and for the meson one M^ ~ 0.790 GeV and ~ 0.14. We see thus that 
the p mass is very close in the large N c limit, to the one in the chiral limit, My. Notice 
that these values for My are far away from those quoted in ref. [13]. The underlying 



reason is that in ref. ][L3|] fy and My were determined directly from the Lagrangian at 
0(p 2 ) in the ENJL expansion, identifying them with their values at q 2 = 0. What we 
find here is that even though the two-point function in eq. Q307] ) has the correct q 2 — > 
limit behaviour it does have, with the choice of vector fields to represent vector particles in 
ref. ||13|| , substantial contributions from higher order terms (mainly of 0(p^) in the ENJL 
expansion). A physical vector field that would include these contributions can in principle 
be defined as is shown by the fact that the inverse of Tly\—q 2 ) is a rather straight line. 
What has happened is that 



^v\-q 2 ) - -— \, : r= ",,^ - (310) 



M 2 (0) - q 2 (l + A + 0(q 2 /A 



The vector meson mass derived in |T3j was My(0) while the slope of the physical two-point 
function (for |g 2 |/A 2 << 1 that is where this ENJL cut-off model makes sense) corresponds 
to rather My ~ My(0) / yl + A. We find from the calculation that indeed A is of order 1 
(A ~ 0.7), explaining the difference in the slope from the 0(p 2 ) ENJL calculation in ref. 
|T3J of the two-point function to the 0(p 4 ) one. 

We can also see from eqs. ( |250| ), ( |262|) - (|264j ) and (|271| ) that the forms of these two- 



point functions are very similar to the corresponding ones in the meson dominance limit 
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but with couplings varying with q 2 . The identification of the corresponding physical values 
will involve analogous procedures to the one described above for the transverse vector 
two-point function one. 



10.3.2 VPP three-point function and the KSRF relation 

In this subsection we discuss how the result for the three-point function nY PP (pi,p2) 



obtained in section |10.1| can be used to determine the physical pion electromagnetic form 
factor in this model. We shall discuss the VPP three-point function flavour structure 
corresponding to the three-point function ITNpi,^) in eq. (|299|) for m = rrii = rrij = rrik 
and p 2 = p\ = p\ for definiteness. 

Since this U+(pi,p 2 ) is a Green's function we first have to reduce the external legs 
to properly normalized pion fields. The vector leg acts here as an external source and is 
properly reduced without bringing in any factor. For this, we first look at the pseudoscalar 
two-point function in eq. ( |264[ ) obtained using the same external fields and parametrize it 
around the pole as 



n, 



-P 



The reducing factor Z n is 



Z w = 



A 2 



(Mi + Mj 



dmU-p 2 ) 



- + hAt (l + 0(m 2 jA 
9s V ~ "V v 



with 



(311) 



dp 



g\{-ml) 



fM - fA-O + K^K, K, o) 



(312) 



and fli 



-q 2 ) in eqs. (g51]) -(g5g) . 



where hipliViil 2 ) defined in appendix and fl(—q 2 
The quantity A is very close to one and exactly one in the chiral limit. Each pion leg brings 
a factor Z\l 2 after reducing the Green's function to the physical amplitude. Rewriting the 
pseudoscalar two-point function in the form in eq. ( |311| ) gives that m 2 is the solution of 
m 2 = m 2 A—m 2 ). 



Reducing the VPP three-point function n+(p 1 ,p 2 ) in eq. (|299|) we find that it can be 
written as follows^ (we shall suppress the flavour indices which are always ii) 



U^( Pl ,p 2 ) 



(p 2 -rrii 



^F VPP (p 2 ,q 2 ) (p 2 -p x y 



(313) 



6 To obtain the j*tt + it three-point function from this is necessary to multiply it by the electric 
charge of the pion. 
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which defines the electromagnetic pion form factor (or in general the pseudoscalar vector 
form factor) Fypp(p 2 ,q 2 ) in this model. (The general pion form factor, i.e different quark 
masses and p\ ^ p\ can be obtained similarly from ( |299|) .) This form factor in the ENJL 
model is expected to be a good approximation at intermediate and low-energy energies, 
within the validity of the ENJL model we are working with, i.e. for \q 2 \ « A 2 ,. The 
explicit expression for this form factorf] is 

1 M 2 (-q 2 



77 / 2 2\ L MV\-1 I U f 2f 2^ 



- i-d - u-mDr-fu-i 2 ) + 2g ? fi 

q — 4m„. 

x \q 2 -2 m l)(jl(-q 2 )~jl(- m l))-4 m %( m lmlq 2 )}} . 

(314) 

Notice that this form factor has no pole at q 2 = 4m 2 . The value of A 2 in eq. (|312|) is 
precisely the one that ensures that Fypp(m 2 , 0) = 1 in the large N c limit as is required by 
the electromagnetic gauge invariance. This must be so since we have imposed the Ward 



identities to obtain this form factor. In figure |Tl] we have plotted the inverse of this form 
factor for the parameters quoted in section |9.9| in the chiral case (m = 0) and in the case 
corresponding to the physical pion mass (m = 3.2 MeV). As can be seen from the picture, 
it is a rather straight line so the complete VMD result for this form factor, i.e., 

F^ P D (ml, q 2 ) = J ^ 7 (315) 

with constant vector mass M p works rather well. The slope of the linear fit of the inverse 
of the form factor in eq. ( |314| ) to this VMD form gives a vector mass which is M p ~ 0.77 



GeV. This mass is very close to the physical value and rather different from the one found 
for the transverse vector two-point function in the VMD limit M p ~ 0.655 GeV in the large 
N c limit. This explains why using the physical p meson mass and the VMD dominance 
works so well but it also shows that this M p "mass" in eq. ( |315| ) has not, in principle, to 
be the same as the mass of the vector meson described by the transverse two-point vector 
function. 

The same three-point function VPP also contains implicitly the p — > tttt coupling 
constant gy. (See Sect. [| for its definition. Notice that is different from the symbol 
defined in ([T6|).) Again, to obtain the physical p — > mr amplitude we should first reduce 
the vector leg that now corresponds to the p particle, (remember that the pion legs have 
been already reduced). This will bring a factor which is similar to the factor 1 + A discussed 
in the previous subsection. We shall, as before, first determine the reducing vector factor 



7 This form factor was also calculated in ref. []57|]. With the appropriate changes of notation it agrees 
with the one found there. 



63 




0.2 1 1 1 1 1 1 1 1 1 1 1 

-0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 

q A 2 (GeV A 2) 



Figure 11: The inverse of the vector form factor of the pion of eq. ( |314| ). For the chi- 
ral limit and with all current quark masses equal to 3.2 MeV. Also plotted is the VMD 
approximation M v (—q 2 )/(M v (—q 2 ) — q 2 ) for the latter case. 



from the vector two-point function in eq. (|307 ). The reducing factor Z p is 



-( 



2f V My 



VyM 2 
B 2 



dM 2 { 



dq 2 



q 2 =Mj, 



(316) 



and is independent 



In this equation the combination 2f v M v is the one given in eq. 
of q 2 . The vector mass M p is again given by the solution to M 2 

One also can rewrite down the electromagnetic pion form factor showing explicitly the 
coupling constant of the p meson to pions, gv, as follows 



M 2 (-M 2 ) 



F V pp — 1 + fvdv 



q 



Mi 



2 ■ 



(317) 



Then, in the complete VMD limit one has fv9v = fllM 2 . In this ENJL model this relation 
is equivalent to gy — (1 — gA)fv-, i- e - one nas complete VMD and the KSRF relation ||111 
2(?y = fv satisfied for = 1/2. 

One can see in the eq. ( |317| ), that reducing the p vector leg brings in a factor B 2 



the numerator and another factor B 2 in the denominator with the net result 



m 

fv(-q 2 )9v(- 
ineq. flHUTp . 



■q 2 ) is not affected by reducing of the vector leg as much as happens to fy{—q 



that 

2^ 
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Figure 12: The generalized KSRF relation. We plot <?y(— g 2 ) for qa = 0.61 (solid line); 
qa — > 1 (short-dashed line) and fv(—q 2 )/2 (dashed line). The difference between the curves 
gives the violation of the KSRF relation. See text for further comments. 



We now define an off-shell coupling gv(—Q 2 ) by Eqs. fl317|) and (|314|) . When expanded 
in q 2 and with m 2 = one gets fy(Q)gy(Q) = 2L 9 , where Lg is the one found in ENJL 
in ref. [I3[, Sect. §. As discussed there, at q 2 = one has the KSRF ||111|| relation, i.e. 



/y(0) = 2gy(0) approximately satisfied. The definition above is the off-shell equivalent to 
the KSRF relation in this model. For qa = the vector mass vanishes and the p meson 
couples as an SU(3)y gauge boson, in fact in this limit one recovers the results of the 
Hidden Gauge Symmetry model [[^J for the non-anomalous sector. In particular, when 
gA = we have that the reducing factor B is 1 as corresponds to external gauge sources. 
In this limit [gA = 0), one still has the KSRF relation analytically satisfied off-shell, i.e. 
f v (-q 2 ) = 2g v (-q 2 ) for all q 2 . 

In the limit gA —> 1 one obtains the constituent quark model result. Let us see how 
gvi—q 2 ) works numerically compared with fv{—Q 2 ) for a definite value of gA- In figure 
12] we plot fv(—q 2 )/2 and gv(—Q 2 ) for the values of parameters discussed in section 



These values correspond to 5u(0) = 0.61. The form factor gv{—Q 2 ) is somewhat dependent 
on q 2 with (2.1 ~ 2.2) gy{—q 2 ) ~ fv{—Q 2 ) m the Euclidean region. In this figure we also 
plot the case gA — ► 1 where the same features can be seen. The form factor gv{—Q 2 ) for any 
value of gA will be between the line fv/2 (i.e., the gA = limit) and the line for gA = 1, 
therefore the KSRF relation is approximately satisfied off-shell for any value of gA- 
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Figure 13: The inverse of the 7r°7*7 form factor for one photon on-shell and one off-shell 
as a function of the photon mass squared, q 2 . Notice the linearity in the Euclidean region. 
Plotted are the full result, M v (-q 2 ) /(M v (-q 2 ) - g 2 )(VMD-like) and the ENJL model 
without vector and axial- vector mesons (qa = !)• 



10.3.3 PVV three-point function 

In this subsection we want to study the 7r°7*7* anomalous form factor. For that we shall 
reduce the PVV Green's function in eq. ( |306| ) calculated in section |10.2| to the physical 



amplitude following the same procedure that in the previous section (for details see there). 
Now, we have to reduce one pion leg, this will bring in a factor y/Z^ and two external 
vector sources legs which are properly reduced without bringing any factor. Then the 
PVV three-point function in eq. ( |306| ) Q can be rewritten as follows 



II- in Nc ,> „V 



x F PVV {q 2 , pIp 2 2 ) (318) 

where Fpyy is the tt° — > 7*7* form factor in this model. Notice that the reducing factor 
A in eq. ( |312| ) goes to one in the chiral limit preserving, in that way, the chiral anomaly 
condition Fpw(Q, 0, 0) = 1. We plot the inverse of this form factor for the case p 2 , = in 
figure |13[ Notice that there Fpyy{m\, 0, 0) 7^ 1 and the difference comes from the reducing 
factor A and is of chiral counting 0(p 6 ). We can expand this form factor for small pf,P2 

8 To obtain the 7r°7*7* three-point function from this II+, is necessary to multiply it by a factor y2 
coming from the ir° flavour structure and a factor e 2 /3 from the quarks electric charge. 
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and pion mass f] as follows 



F PVV (ml,pl,pl) = l + p{pl + pl) +p 'ml + 0(q 4 ) 
this expansion defines the slopes p and p' which in this model are 

1 r(2,M 2 /A 2 )\ 



(319) 



and 



P 



P 



9a(0) 



9a(0) 



M^(0) 12M 2 J ' 

T(2,M 2 /A 2 ) r(l,M 2 /A 2 .; 



12M 2 



12r(0,M 2 /A 2 )M : 



(320) 



Where the second term in p' comes from the reducing factor A defined. The constituent 
quark mass M here is the one corresponding to the current quark mass value m = 3.2 MeV 
used in the numerical applications section |9~9| . Using My-(O) = 6M 2 g J 4(0)/(l — 5u(0)) Jl3| 
we can write down them as 



and 



P 



P 




2-T(2,M 2 /A 2 ) WO) 



r(2,M 2 /A 2 ; 



r(i,M 2 /A 2 )' 
r(o, m 2 /a 2 ) / 



(321) 



which interpolate between the constituent quark-model result #^(0) = 1 and the gauge 
vector meson result <?a(0) = 0. 

With the input parameters we have been using (see numerical application section |9.9| ) 
we get 



and 



p = (0.86 + 0.67) = 1.53 GeV" 2 
y = (0.67 - 0.27) = 0.40 GeV" 2 . 



(322) 



Where for p the first number between brackets is the vector meson exchange contribution 
and the second is the constituent quark contribution (up to gr^(0)). We see that both con- 
tributions are very similar giving some kind of complementarity between both approaches 
and explaining the relative success of both when used to describe this slope. For p' they are 
the constituent quark contribution and the one coming from the pion leg reducing factor 
1/A (Notice the cancellation there.) Experimentally [[L12|| 



p = (1.8 ±0.14) GeV" 



(323) 



Taking into account that the 1/N C corrections from x?T loops are estimated [ 1 1 3 1 to be 
twice the experimental error we consider the result as good. The cancellation in p' is also 



3 For the ir° decay we are on the pole and hence q 2 
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welcome since otherwise the SU(3) breaking corrections in the rj decay would have been 
much too large. 

Let us compare this full result in eq. (|321| ) with the one obtained in ref. [61] in this 
same model assuming complete VMD in the chiral limit. There, the same prescription to 
include the QCD chiral anomaly that here was used at the one-loop level with the 
result 



P = T^^P- (324) 



12M 2 g A (0) 

Of course, this complete VMD result vanishes when qa = 1 where vector mesons decouple. 



11 Hadronic Matrix elements 



11.1 The 7T+ - tt° Mass Difference 



In this section we will discuss the general philosophy behind the 1/N C method of calculating 
nonleptonic matrix elements. A good review where also the references to the original 
papers can be found are the lecturers by Gerard fll4||. The application of this method to 



the 7r + — 7T° mass difference can be found in Ref. [ |115|1 and the calculation within the QCD 
effective action model and the EN JL model is in Refs. [0, [53] . 

We look at this quantity because it is the simplest nonleptonic matrix elements in 
several respects. There is no factorizable contribution because the photon is spin 1 and the 
pion spin 0. It involves only pions so we expect the limit where the current quark masses 
vanish to be a good approximation and (unlike Br) it doesn't vanish and is well defined 
in this limit. The latter remark has one very useful consequence. Using PCAC it can be 



shown [|116|| that this matrix element can be related to a vacuum matrix element. So the 
mass difference becomes a vacuum matrix element of the photon propagator integrated 
over all momenta in the presence of the strong interactions. Schematically, the matrix 
element (tt + \ J 2 \tt + ) can be rewritten in terms of (0| J 2 |0). The precise expression in terms 
of the hadronic two-point functions is given by [|116| , |115| , |7] 



2 2 3ft em 



fdQ*Q*(nPm-rt}\QZ)) . (325) 



Svr/ 2 

Eq. (|325|) involves an integral over all distance scales. The underlying idea is now to split 



this integral into two parts, = Jo + 1^2, and then to evaluate both pieces separately. 

The long distance part in 1/N C can be calculated in models since in 1/N C the only quan- 
tities needed are the couplings of currents to hadrons and not of full four-quark operators 
to hadrons. The essence of the 1/N C method is to do the short-distance part using the op- 
erator expansion and then use 1/N C to evaluate the matrix element. Here this corresponds 
to using as the difference of 2-point functions: 

(nf(Q 2 ) - n?(Q 2 )) = =±8na s (Q>) {{qq)f . (326) 
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This leads to[ |115| 



Ami 



SD 



3aas 



my 



(327) 



One can then still do a renormalization group improvement of this|53] . 

The long distance part of the integral requires more care. There are several approaches. 

1. One can take the measured spectral functions and use these to evaluate the two-point 
functions needed in the integral. The most recent evaluation of this is in Ref. ||117 . 



2. The two-point functions can be approximated by including the p, tt and ai contri- 



bution. This was done neglecting the QCD part in the original paper ||116|1 and more 
recently in ||115| . 



3. We can take only the 7r contribution ||1 15| . This is most like the original 1/N C 
method(Ref. ||114j| and references therein). This leads to AmJ = 3a/i 2 /(47r). 



4. One can use the QCD effective action approach] 

5. The ENJL model can be used I 



All of these approaches give a good result for the mass difference. In cases 1,2 and 5 a 
good matching was also obtained. This means that we can vary /i, the split between the 
short- and long-distance part of the integral, over a reasonable interval without changing 
the result. In figure [14] the long-distance result with only the pion is shown and the ENJL 
long-distance result. Also shown is the experimental value, the short-distance result and 
the sum of short- and long-distance for the ENJL case. The value of (qq) used is the one 
given by the ENJL model. 

At this point I would like to remark that for this quantity in the QCD effective action 
approach one only obtains a gauge invariant result if the pion is explicitly taken as propa- 
gating (see |7[]). This shows that in this model the pion degree of freedom has to be added 
by hand. The gauge dependence then cancels between a two- and a three-loop diagram. 



11.2 B 



K 



In this section the extension to weak nonleptonic matrix elements of the methods in the 
previous section is discussed on the example of Bk- Here again the pure 1/N C method |114| , 
118 , the QCD effective action model|50| and the ENJL model|59], [60| . An overview of 
theoretical situation a few years ago can be found in Ref.[|119||. The main alternatives to 



the present method are lattice calculations [12C, E| and 2 and 3-point QCD sum rules[121 



The short-distance integration here is done using the renormalization group. This sums 
the possible large logarithms involving log(m^//i 2 ). The problem then reduces to the study 
of 

(K°\O A s = 2(x)\K ) = -B K (ji)J%m a K (328) 
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Figure 14: The results for m 2 r+ — m 2 ^: long-distance result with only the pion (LD- 
CHPT); ENJL long-distance (LD-ENJL); experimental value (exp.); the short- distance 
result (SD281) and the sum of short- and long-distance ENJL (full). 



with the AS = 2 operator AS=2 (x) = Lf(x)L^ sd (x); 2Lf(x) = ^(x)^ (1 - 75) d(x) and 
summation over colours is understood. Eq. ( p28| ) is also the definition of the B^ parameter. 
The different approximations give 

1. Vacuum Insertion : Bk((i) = 1. 

2. Leading in 1/N C : B K {fi) = 3/4. 

3. The standard 1/N C result: B K (p) = 3/4(1 - 2/i 2 /(167r 2 /*))- 

4. 1/N C with inclusion of vector mesons |il4|| : 



2 ,,2 



7 2 3 my/I 



4 \i 2 + m 



3 2 // 2 + mf/ 
-m v log ■ 



mi 



(329) 



5. The QCD effective action result [pO 

3 



(—G 2 ) + ... 



(330) 



We would also like to study the effects of off-shellness. Therefore we do not directly study 
the matrix element in Eq. ( |328|) but the Green function 



Gf nAs= 



d 4 xe iq - x (0\T (P ds (0)P ds (x)T AS=2 ) |0) 



(331) 
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in the presence of strong interactions. We use the ENJL model for scales below or 
around the spontaneous symmetry breaking scale. Here Gp is the Fermi coupling con- 
stant, we use P ds (x) = d(x)ij 5 s(x), with summation over colour understood and Tas=2 = 
— Gp J d A y 0/\s=2(y)- The reason to calculate this two-point function rather than directly 
the matrix element is that we can now perform the calculation fully in the Euclidean region 
so we do not have the problem of imaginary scalar products. This also allows us in princi- 
ple to obtain an estimate of off-shell effects in the matrix elements. This will be important 
in later work to assess the uncertainty when trying to extrapolate from K —>■ n decays to 
K — > 27T. This quantity is also very similar to what is used in the lattice and QCD sum 
rule calculations of B K . 

The AS = 2 operator can be rewritten as 

T AS =2 = -G F / d 4 x, J d'x 2 e-^-^L^ Xl )L^x 2 ). (332) 

This allows us to consider this operator as being produced at the Mw scale by the exchange 
of a heavy X AS = 2 boson. We will work in the Euclidean domain where all momenta 
squared are negative. The integral in the modulus of the momentum r in ( p32j ) is then split 
into two parts, J^ w d\r\ = $ d\r\ + f™ w d\r\. In principle one should then evaluate both 
parts separately as was done for the tt + — ir° mass difference in the above quoted references. 
Here we will do the upper part of the integral using the renormalization group. This results 
in the integral being of the same form but multiplied with the Wilson coefficient C(n), 

T AS =2 = —Gp C(ji) jT -0- J d 4 Xl J d 4 x 2 e- ir ^-^L s «( Xl )L>: d (x 2 ) . (333) 

This can now be studied using the 1/N C expansion. We can first do this within a chiral 



expansion leading to the result in the chiral limit (see Ref.|59| for details): 

(334) 



3 / 1 

Bk(h)chpt — 7 1 



4 \* 16tt 2 / 2 



^4 



The correction is negative. It disagrees somewhat with the result obtained in [ |118|| because 
there no attempt at identifying the cut-off across different diagrams was made. Since we 
work at leading level in 1/N C in the NLO CHPT corrections we have included the relevant 
singlet (?7i) component as well using nonet symmetry. The correction in ( |334j ) has precisely 
the right behaviour to cancel partly C(fi) which increases with increasing /i. 

The same calculation can now be performed for the ENJL model. Here the major 
complication is the number of different diagrams that has to be evaluated. An example of 
one of the classes is shown in figure |Tj| 

We now evaluate all contributions numerically to the two-point function of Eq. ( |33 1| ) . 
The results for several input values are in table |[ We have studied three cases, namely, 
the chiral case, rrid = m s = 0, the case with SU(3) symmetry breaking m s = 83 MeV ^ 
rrid = 3 MeV and the case with m s = rrid = 43 MeV. The other parameters are G,s = 
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Figure 15: A leading 1/N C contribution to the nonfactorizable part of Has=2{<1 2 ) in the NJL 
model. The crosshatched areas are the full two-point functions as discussed in subsection 
3.1. Point E and F are connected via T&s=2- 



1.216, A x = 1.16 GeV and Gy = 0. The latter simplifies the calculation by about an 
order of magnitude. Preliminary results for the Gy ^ case have the same qualitative 
conclusion |K| but typically somewhat lower values of Bk and less good matching. 

The procedure we have followed to analyze the numerical results is the following. We 
fit the ratio between the correction and the leading 1/N C result for a fixed scale /i to 
a/q 2 + b + cq 2 which always gives a very good fit (a, b and c are // dependent). Once we 
have this fit we can extrapolate our Bk form factor (remember that we have calculated it 
for Euclidean q 2 ) to the physical Bk, i.e. to q 2 = m 2 K ~ 0,0.13 GeV 2 (chiral,other cases). 

Let us first treat the chiral or massless quarks case. Here a nontrivial check on the 
results is that the diagrams have a behaviour which sums to 1/q 2 , i.e. a should be zero. The 



Table 2: Results for Bk and Bk in the ENJL model. 



fjt (GeV) 


b x M 


n K 




B a M 


DtU 

°K 




D K 


0.3 


0.68 


0.50 


0.74 


0.50 


0.55 


0.74 


0.55 


0.5 


0.59 


0.59 


0.71 


-0.44 


0.71 


0.72 


0.72 


0.7 


0.53 


0.58 


0.69 


-2 


0.75 


0.68 


0.75 


0.9 


0.48 


0.55 


0.66 


-3 


0.76 


0.65 


0.75 


1.1 


0.45 


0.54 


0.64 


-4 


0.76 


0.64 


0.76 
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individual contributions do not have this behaviour, b is the relevant contribution to Bk 
since m 2 K \ x = 0. The first three columns in table || are //, B K (n) and B\ = B K (fi)as(f^) a+ 
with a + = —2/9 and = 250 MeV. The hatted quantity is the scale independent 

quantity. Good matching is obtained if this value is stable within a range of \i. 

In the second case (m s ^ rrid), because the chiral symmetry is broken, there is a 
possibility for contributions to Bx(q 2 ) that are not proportional to q 2 , i.e. a / 0. In fact 
a CHPT calculation predicts precisely the presence of this type of terms pDf. For small 
values of q 2 the part due to a dominates even though it is only a small correction when 
extrapolating to the physical B^ at q 2 = m 2 K . This can be found in column 4. The fifth 
column is the form factor B K for q 2 = —0.001 GeV 2 where the correction due to the a term 
is sizeable. Notice the difference between these two columns. This same feature should 
be visible in the lattice calculations as soon as they are done with different quark masses. 
The invariant B^ for this case is in column 6. 

In the the last case, i.e. = m s , which is similar to the present lattice QCD calcu- 
lations, the fit gives a = to a good precision and the value of B e £ extracted is rather 
independent of q 2 . The invariant B e ^ in this case is in column 8. 

In view of the results of O, ITDL |54| we expect to get a good prediction for the effects 



of non-zero and different current quark masses. We see those and find a significant change 
due to both: 

B%(m 2 K « 0.13 GeV 2 ) « 1.35 B K {m 2 K = 0) (335) 

for scales \x ~ (0.7 ~ 1.1) GeV. For the extrapolation to the kaon pole the difference 
between the masses has a much smaller effect than the fact that they were non-zero. In 
order to compute Bk in the general case a careful extrapolation to the poles was needed. 
The final correction to the Bk parameter compared to its leading value of 3/4 turns out 
to be rather small. 



12 Conclusions 

In this Physics Reports an overview of the results and the methods of references 



51, pq, B4l K9j was given. There have also been several short versions given in talks and 



lectures JB6|, |65], 0. The main conclusion is that it is possible to use constituent chiral 
quark models and in particular the extended Nambu-Jona-Lasinio models to obtain results 
in the hadronic sector. The method underlying most of the results reviewed here allows 
for a clean separation of cut-off effects and consequences of the general structure of the 
model. This was especially clear in the case of the low-energy expansion, section || and 
the two-point functions, section [|. It also allowed for a discussion of the anomaly in the 
presence of extra pointlike quark interactions, section [j]. 

In the case of 3 and higher point functions the number of free functions in the general 
approach rapidly grows but still some general features are visible by the reduction to the 
basic one-loop functions. 

The results also allowed the derivation of relations that are valid in a large class of 
Nambu-Jona-Lasinio like models. In particular the one with the inclusion of gluonic cor- 
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rections. In fact the final results are rather insensitive to the inclusion of these extra effects 
precisely because the success of the ENJL model rests to a large extent on the relations 
we have derived here. These are also well satisfied by the hadronic experimental quantities 
we have tried to explain with this approach. 

Especially the presence of several of the short-distance constraints from QCD in the 
ENJL model makes it a prime candidate for trying to estimate semi-analytically hadronic 
matrix elements. This has been done at present for the tt + — 7r° mass difference and the 
Bpc parameter. Both results have been reviewed here shortly. 

In conclusion the simple extended nambu-Jona-Lasinio model provides us with a com- 
pact way to describe low-energy hadronic physics and understand a large body of ex- 
perimental results. It also provides a framework to systematically explore effects due to 
explicit chiral symmetry breaking beyond chiral perturbation theory as shown here in the 
discussion about the Bx parameter. 
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A Derivation of the Ward identities 



In this appendix we generalize the proof in the appendix of ref. to the case with 
nonzero current quark masses. There a proof was given of all relevant identities in terms 
of the heat kernel expansion (for an excellent recent review and definitions see ref. ||45|| ) 
and some of them in terms of the Ward identities as well. Here those which can be derived 
directly from the Ward identities can also be derived from the heat kernel expansion but 
since they involve different masses they require a resummation of different terms. For these 
the direct derivation of the Ward identities is actually simpler. Only for the additional 
relations will we give the heat kernel derivation. 

At the one-loop level we use as Lagrangian the one in eq. flHf ) (with the same definitions 
as there) 

C ENJL = qDq (A.l) 

where V contains the couplings to the external fields Z^r^s and p as well as the effects 
of the four-quark terms in /^'jl anc ^ ^njl on the Q uar k currents at the one-loop level. 
In particular it contains the constituent quark masses, Mj. However, we shall keep the 
notation l^r^s and p to denote the quark current sources in the presence of these four- 
quark NJL operators. The one-loop current identities derived from this Lagrangian are 

WV* = -i (M< - M,-) S* j 

= (Mi + Mj) P ij . (A. 2) 



74 



When the whole series of constituent quark bubbles are summed these identities are sat- 
isfied changing constituent quark masses by current quark masses. In addition we use the 
equal time commutation relations for fermions 

{(#(*), q> p (y)} _ = *<WM 3 (* - y) • (A.3) 

Here a and /3 are Dirac indices and x means the spatial components of x. Multiplying the 
two-point functions with iq^ is equivalent to taking a derivative of the exponential under 
the integrals in eqs. (|213|) to ( |216|) . By partial integration we then get several terms, 



those due to the time ordering which leads to equal time commutators and those where 
the derivative hits one of the currents. The first type are evaluated using eq. ( |A.3|) and 



the second type are related to other two-point functions using eq. (|A.2| ). This then leads 
to the expressions ( |234| ) to ( |235| ). 

The derivation of the other two identities is slightly more complicated. The effective 
action of the Lagrangian in eq. ( |A.1| ) can be obtained in Euclidean space as a heat kernel 



expansion (see ref. f45[). The coefficients of this expansion are the so-called Seeley-DeWitt 
coefficients, they are constructed out of the two quantities E and R^. These are defined 

as 

V ] V = - V M + E + M 2 , 
= [V^V,] , 

V M # ee ^#-i[i; M ,#]-i[a M 75,#]. (A.4) 
If in eq. (|A. 1| ) the Dirac operator T> contains couplings to gluons these should not be 



taken into account in eq. ( |A.4j ). The relevant heat kernel expansion in that case will have 
different coefficients depending on vacuum expectation values of gluonic operators, but will 
still be constructed out of the quantities in eq. ( |A.4j ) (depending now also on the gluon 
field) . The quantity M is the mass that is used in the heat kernel expansion. The operator 
V is 

iY{d^ - iVf, - ia^ 5 ) - M - s + ip^ 5 . (A. 5) 

Here Ai = diag(m u , rrid, m s ) is the current quark mass matrix and we allow for spontaneous 
chiral symmetry breaking solution (0|s(x)|0) 7^ 0. For the terms relevant to two-point 
functions we have 



and 



E = -'-a^R^ + tYd^M + s + tm.) 

- Y {a M 7 5 , M + s — ip 75 } + {M,s}-i [M, p] 7s + s 2 + p 2 

+ M 2 — ~M 2 , 

with 

tf ee av-av-i[/,/], 

oT ee & k a v - d v a" - i[a», a v ] , 

d^# ee c^#-zK,#]. (A.6) 
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The main difference with ref. is the occurrence of the last line in the expression for 
E in (|A.6|) . We shall call this last line E . In this equation, M = diag(M u , M d , M s ), 
the diagonal matrix of the constituent quark masses defined in eq. fl20|). Notice that the 
scalar field here has been shifted and we have now (0|s(x)|0) = (though we use the 
same notation for it). When G,s — > in eq. ( pOp then M — > and M — > A4. Let us 
now systematically go through all possible types of terms in the expansion. We shall not 
discuss the mixed two-point functions here since we only want to prove eqs. ( |238| )-( p39| ). 
In the heat kernel expansion, those terms containing two factors only contribute 

to the transverse parts, Uy A and in the same way. Their contributions hence obviously 
satisfy eqs. ( E38 )-( E3^ ). Similarly, one factor R^ u requires the presence of two covariant 
derivatives V M . By commuting derivatives (the extra terms only contribute to three and 
higher point functions) and partial integration these can be brought next to each other so 
they convert into a second factor R^ u . This brings us back to the previous case. Intervening 
_E"s can only contribute via E but these do not spoil the above argument. The first term 
in E, namely a^R^, requires a 2nd a^R^ because otherwise the trace over Dirac indices 
vanishes. These also behave like terms with two factors R^ v . Therefore, in the remainder 
we are only concerned with E without this first term. 

E can also directly contribute to the scalar and pseudoscalar two-point function in the 
same way via s 2 + p 2 . Extra factors E become again Eq and extra derivatives also respect 
the relation ( |239|) . The most complicated case is where both fields come from a different 
E. This contributes in the form EqEE^O^E. These contribute to all form factors in the 
form M™ M™ q 21 times the coefficients listed in Table |3|. These coefficients obviously satisfy 



Function 


Contribution 


rip 

TT (0) 

TT (1) 
n (o) 

n (1) 


-q 2 + (Mi + Mj) 2 
-q 2 + (Mj - Mj) 2 

(M. + Mrf/q 2 
-{Mi + Mj) 2 jq 2 

(Mt-Mrf/q 2 
-( Ml -M 3 ) 2 lq 2 



Table 3: The contribution of terms of the type E m+n+2 to the two-point functions. 

the relations ( P38| )- ([239|) . The last type of terms is where the external fields come out of a 
derivative. We do not consider the mixed case here, so both the fields have to come out of 
a derivative due to the 7 M that is necessarily present in the E that would be a candidate for 
the external field. So there are those where the external fields are contained in two factors 
V M . If the indices of these are different, then there need to be at least two extra derivatives 

present that will produce a q^qv This contributes equally to Uy^ and n^. If the indices 
are equal, it will contribute proportional to g^ v and thus to the vector and axial-vector 



76 



equally with = 0. This completes the proof of the identities ( |238| )-([23 

Now it remains to prove that these contributions will never produce a pole in n ( ' 0+1 ' ) at 
q 2 = 0. Terms that contain two factors R^ v contain two factors of momenta and hence do 
not. Terms with one factor R^ v can be brought in the form with two so do not produce 

a pole either. From Table |^ there is no contribution from that type of terms to Ily^ 1 . 
Then those with external fields from with different derivatives necessarily contain extra 
factors q^q v so do not contribute to a possible pole at q 2 = and the last type of terms 

does not contribute to HyA as shown above. This completes the proof. 

B Explicit expressions for the barred two-point func- 
tions 

Here we shall give the one-constituent-quark-loop expression for the two-point functions 
defined in eqs. ( 213 )- ([218| ) in the presence of current quark masses. These two-point 



functions are denoted in the text as the II ones. They fulfil the same Ward identities as 
the full-ones in eqs. (|234|) -( |23~7D changing the current quark masses there by the constituent 
quark ones. In addition, they also satisfy the Ward identities in eqs. (|238|)-(|2"3Tf). Using 



these identities one can see that there are only two independent functions out of ITy ^ , ILy , 

11^ , , Il!!f Tip , lis and lip. We shall take lip and Tly^ + Uy as these functions. The 
explicit expressions are 



(n£ } + 4 0) ) {Q% = —^8j o dxx(l-x)T(0,x tJ ), (B.l) 
Kp(Q 2 k = T^ £ 2 4/ 1 dx(M 4 x + M J (l-x))r(0,x u ), (B.2) 

l07T JO 



where 



M 2 x + M 2 (l-x) + Q 2 x(l-x) 

Xij — ^ 2 . l-D-<jJ 

X 

One can obtain all the others one-loop two-point functions in function of these two by 
using the Ward identities mentioned above. For instance, for the lly one gets 



(o) 2 (M.-MjY Up (Q )jj 

] U M i + M J Q 2 (Q 2 + (M i -M j ) 2 ) (B.4) 



x {(Mi + M 3 f + 9A (Q*hrnUQ>) (l - (w^) + Q 2 } 
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C Explicit expression for the one-loop form factor 

n^(pi,p 2 ) 

Here we shall give the one-constituent-quark-loop expression for the three-point function 
U. (pi,P2) defined in eq. ( |285| ). We shall give it for Mi = Mk = M m . The explicit 
expression is (remember that we have j = m), 



+ 



1 (n?(-P?)« 
Pl ' P2 ^(np^A-npW 



2M, 



PlP2 ~ (Pi • P2 

+ (pi • P2) (n p (-Pi)u - Hp (-g 



2 

+ (Pi <-> 



— h(pi,pl,q 2 )pl 



-P2) • 



, . / \ Pl + (Pl-P2) 
1 + (Pi • P 2 ) — 



PlP2 - (Pi ■ P2) . 



(C.l) 

Where the two-point function 11^ (— p 2 ) was given in appendix [B| and the function Ij,{p\-,p\i <Z 2 ) 
is 



N r 1 r 1 

h(plpl,q 2 ) = —^2M? dxx dy 

l07T JO JO 



rl /"I r(l,M 2 (x ?2 /)/A 2 ) 



M 2 (x,?/) 



(C.2) 



with 



M 2 (x,y) 



M/ - p 2 (l - x) - p\x{\ -y) + (pi(l - x) - p 2 z(l - j/))' 



(C.3) 
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